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CSl ■ ABSTRACT 

Aims. We have performed a comprehensive parameter study of the collapse of rotating, strongly magnetized stellar cores in axisymmetry to 
determine their gravitational wave signature based on the Einstein quadrupole formula. 

Methods. We use a Newtonian explicit magnetohydrodynamic Eulerian code based on the relaxing-TVD method for the solution of the ideal 
MHD equations, and apply the constraint-transport method to guarantee a divergence-free evolution of the magnetic field. We neglect effects 
due to neutrino transport and employ a simplified equation of state. The initial models are polytropes in rotational equilibrium with a prescribed 
degree of differential rotation and rotational energy. The initial magnetic fields are purely poloidal the field strength ranging from lO'^G to 
' lO'^G. The evolution of the core is followed until a few ten milliseconds past core bounce. 
QQ Results. The initial magnetic fields are amplified mainly by the differential rotation of the core giving rise to a strong toroidal field component 
y—( with an energy comparable to the rotational energy. The poloidal field component grows by compression during collapse, but does not change 
significantly after core bounce. In large parts of the simulated cores the growth time of the magneto-rotational instability (MRI) is of the 
T-H ^ order of a few milliseconds. The saturation field strengths that can be reached both via a pure CI dynamo or the MRI are of the order of 
. lO'^G at the surface of the core. Sheet-like circulation flows which produce a strong poloidal field component transporting angular momentum 
■ outwards develop due to MRI, provided the initial field is not too weak. Weak initial magnetic fields 10" G) have no significant effect on 
^ the dynamics of the core and the gravitational wave signal. Strong initial fields (S lO'^G) cause considerable angular momentum transport 
O I whereby rotational energy is extracted from the collapsed core which loses centrifugal support and enters a phase of secular contraction. The 
gravitational wave amplitude at bounce changes by up to a few ten percent compared to the corresponding non-magnetic model. If the angular 
momentum losses are large, the post-bounce model. If the angular momentum losses are large the post-bounce equilibrium state of the core 
changes from a centrifugally to a pressure supported one. This transition imprints in the gravitational wave signal a reduction of the amplitude 
^ of the large-scale oscillations characteristic of cores bouncing due to centrifugal forces. 

• • 1 In some models the quasi-periodic large-scale oscillations are replaced by higher frequency irregular oscillations. This pattern defines a new 
^ ■ signal type which we call a type IV gravitational wave signal. Collimated bipolar outflows give rise to a unique feature that may allow their 
detection by means of gravitational wave astronomy: a large positive quadrupole wave amplitude of similar size as that of the bounce signal. 
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1. Introduction 

The gravitational binding energy liberated by the collapse of 
the iron core of a massive (M ^ 8 Mq) star to a neutron star 
is the commonly accepted energy source of type Ib/c and type 
II supernovae, as a few percent of this energy are sufficient to 
unbind and rapidly eject the stellar envelope and to create the 
supernova outburst. However, which physical processes turn 
the central implosion into the explosion of the stellar layers 
suiTounding the forming neutron star is still debated in spite 
of many efforts over more than three decades. Heating of stel- 
lar gas just outside the proto-neutron star (PNS) by neutrinos 
diffusing and being advected out of its interior is thought to 
play a crucial role in the explosion mechanism. However, as 



current neutrino-driven supernova models produce (weak) ex- 
plosions only for low mass progenitors (for a recent review, 
see e.g. Janka et al. 2004), there may be a need to include addi- 
tional physics in the models in order to make them work suc- 
cessfully for more massive progenitors, too. 

On this account magneto-rotational core collapse, 

which has been s t udied b y a few authors in t h e past 
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iMeier et al . 1976; Muller & Hillebrandt 1979; Ohnishi 1983; 
Svmbalist vi 19841) . has become an active research field in recent 
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the necessarily multi-dimensional magneto-hydrodynamic 
(MHD) sim ulations, observations indi c ating very asynimetri c 
explosions dWang et all Il996l l200lt Leonard et all 1200 iL 
and the interpretation of Anomalous X-Ray Pulsars and 
Soft Gamma-Ray Repeaters as magnetars, i.e. very s t rongly 

magnetized neutron star s j Duncan & ThomDson| Il992t 

iThompson & Duncanll996tlKouveliotou et al...l999t) . 

Concerning the initial conditions for magneto-rotational 
core collapse the up to now most ad vanced ev olutionary cal- 
culations of rotating massive stars faeger et al.. .2005ii) pre- 
dict that the initial rotation rates are more than an order of 
magnitude smaller than (i) the minimum ones used in past 
(parameter) studies of magneto-rotational core collapse, and 
(ii) those predicted by previous evolutionar y calculations (see, 
e.g. lWooslev et al J2OO2;' Hi rschi et al J2OO3) which lead to neu- 
tron stars rotating very rapidly (~ 1 ms) at birth. The latter stud- 
ies ignored the torques exerted in differentially rotating regions 
by the magnetic fields that thread them. Thus, the stars end up 
with 3 to 50 time s more angular momentum than in the mod- 
els by iHeger etal] ll2005h in that part of their core destined to 
collapse to a neutron star. 

The strength (and distribution) of the initial magnetic field 
in the stellar core is unknown. If weak initially, several possi- 
ble amplification mechanisms exist that may amplify the mag- 
netic field of the collapsing progenitor to a dynamically im- 
portant strength. Linear amplifica tion of the field by means 
of diff'erential rotation will occur iMeier etaP ( Il976l) . which 
transforms rotational energy into magnetic energy by wind- 
ing up any seed polodial field into a toroidal magnetic field. 
This process can be accompanied by the action of meridional 
(e.g. convective) motions that transform toroidal into poloidal 
fields. Both processes together lead to the so-called a-Q. dy- 
namo. R ecently, th e magneto-rotational instability (MRI) (see 
iBalbus & Hawleylil998) has received a lot of interest in the 
context of supernova collapse and explosion (lAkivama et al 



2003'; Kotake et al. 2004a: Vamada & Saw ai"2004': Sawa i et al 



2005) . Unlike linear wrapping, the MRI will give rise to an 
exponential growth of the fie ld strength while w orking on the 
same time scale (see however lSawai e t al.''2005V The MRI sat- 
uration field is independent of the initial field, i.e. even quite 
small initial fields can be amplified to dynamically important 
strengths. The MRI will occur if the radial gradient of the an- 
gular velocity is negative, a condition arising quite naturally in 
core collapse situations. 

A major effect of magnetic fields on the collapse dynamics 
is the transport of angular mome ntum. Due to its v ery low (fluid 
and V shear) viscosity (see e.g. iKeil etal.llT99fil) a collapsing 
non-magnetized stellar core maintains its Lagrangian angular 
momentum profile j(m), m being the Lagrangian mass coordi- 
nate, on time scales of ~ Is, but magneti c fields can signifi - 
cantly redistribute the angular momentum (Meier et al. '1976^. 
This can slow down the forming neutron star and thus counter- 
act the eff'ects of rotation. In some cases, even retrograde rota- 
tion may result in some parts of the core ("Muller & Hillebrandt' 
Angular momentum transfer can also destabilize the 
rotational equilibrium the core resides in after a centrifugal 
bounce at sub-nuclear densities, and lead to a subsequent (sec- 
ond) collapse to nuclear densities and beyond that releases 



large amounts of gravitational binding energy dSvmbalistvl 
Il984l) . Conversely, the violent convective flow both inside the 
neutrino sphere and between the neutrino sphere and the shock 
will transport and amplify magnetic fields in the collapsed core 



Analytic considerations (Meier etal. 1976; 
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i2005i) show that magneto-rotational core collapse might lead 
to jet-like explosions. Though the magnetic stress will remain 
below equipartition strength in most regions of the star, it might 
aff'ect the dynamics of the core through its anisotropic compo- 
nents. Magnetic stresses can assist in pushing the stalled shock, 
or may even drive a mildly relativistic outflow in form of a jet 
along the rotation axis, which is powered by the rotational en- 
ergy transfered to the jet by the magnetic stresses. 

Observations of gravitational waves (GWs) will al- 
low one to learn more about the supernova mecha- 
nism as they provide pristine information directly from 
the stellar interior, in particular about the amount and 
distribution of the angular momentum, and the strength 
and topo l ogy o f the magnetic field . Simp l ified Newtonian 



(Muefler 1982 
. 1994: 2: werger & Muller. 
iFrver et al .ll2(XMr [ 



Monchmeyer et al 
1997. . 

IBJ\ 



^'"f etal. I 



1991 
Kotake et al 
[Yamada & 



Yamada & Sato* 




general-relativistic (Di mmelmeier et al . 2002a b) calculations 
of rotational core collapse predict the emission of a strong sig- 
nal around core bounce, and that the magnitude of the bounce 
signal as well as the post-bounce gravitational radiation de- 
pend sensitively on the initial rotation rate and rotation pro- 
file. Newtonian hydrodynamic simulations using more sophis- 
ticated micro- and tran sport-physics, as w ell as state-of-the- 
art rotating progenitors dHeger et alJ2005l) show that the gravi- 
tational wave signal at core bounce is small compared to the 
signal produced by convective motions in the post-bounc e 
core and by aspheric neutrino emission iMiiller et alJ ll2004l) . 
Magneto-rotational effects on the g ravitational w a ve sign ature 
were first inves tigated in detail by Kotake et al. (2004b) and 
I Yamada & Sawai (,2004) who found differences from the sig- 
nature of purely hydrodynamic models only in the case of very 
strong initial fields (\B\ > lO'^G). 

In the following we present a comprehensive parameter 
study of the axisymmetric Newtonian core collapse of ro- 
tating magnetized polytropes and of their gravitatio nal wave 
signatu re. Our study extends the work of Zwerger & Miilleil 
(I1997I) (hereafter ZM) and the complementary work of 
'Dimmelm eier et aD ('2002a"b) (hereafter DFM), who investi- 
gated the hydrodynamic collapse of a large set of rotating poly- 
tropes in Newtonian and general relativistic gravity, respec- 
tively. The simulations have been performed with the recently 
developed MHD difference scheme of Pen et al. (2003). They 
incorporate neither neutrino transport nor nuclear burning pro- 
cesses. Due to the reduced complexity of our models we could 
explore many of them covering a large region in parameter 
space, the focus being the gravitational wave signal emitted by 
magnetized stellar cores and their dynamic evolution. Because 
of our assumptions and approximations the vahdity of our mod- 
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els is limited to the stage of core collapse and to the first few 
ten milliseconds of their post-bounce evolution. 

Several related but less comprehensive numerical MHD 
studies have be en performed in the past few years: 
lYamada & Sawai (2004) used th e ZEUS-2D code, emp loyed 
the parametric equation of state o f lYamada & Satollligg-A . con- 
sidered no neutrinos, and followed the evolution of initially 
rapidly rotating and very strongly magnetized ^ lO'^G 
cores with a purel y homogeneous poloidal field. Contrary to 
iLeBlanc & Wilsoni(d97a) and .S vmbaUstv. (J 984) they find that 
the magnetic field becomes strongest behind the shock wave 
and not in the inner core, and thus is the main driving factor 
of the observed jet outflow along the rotation axis. Besides a 
field amplification by differential rotation, they also observe 
the possible action of the MRI. They calculate the gravita- 
tional wave signal in the quadrupole approximation finding 
no substantial difference between t he bounce signal of m agne- 
tized and non-magnetized models. FKotake etall (|2004b) also 
use the ZEUS-2D code to which they add an approximate 
neutrino cooling with a leakage scheme. They assume an ini- 
tially predominantly toroidal magnetic field in their investi- 
gated 14 models of which all but one are very strongly mag- 
netized \B\ > 10 ^ 'G. Be sides the simplified equation of state of 
lYamada & Sato (1994') they also co nsider two realistic equa- 
tions of state. Kotake et alJ ll2004bl) focus their study on the 
effect of the magnetic field on the gravitational wave signal, 
and find that the gravitational wave amplitudes are lowered by 
~ 10% for mo dels with the stronge st initial magneti c fields 
(Ifll ~ lO'Vi-) iKotake et alJ ll2004ah . iTakiwaki et aP (|2004), 
and Kotake et al. ' ('2005') all using the same input physics and 
numerics as Kotake et al, ( „2004b) are concerned with the ef- 
fects of the magnetic fields on the anisotropic neutrino radia- 
tion and convection, on the propagation of the shock wave, and 
on the rotation-induced anisotropic neutrino hea t ing thr ough 
parity-violating effects, respectivelv. lKotake et al ] |l2004ah find 
that the aspherical shapes of the shock and of the neutrino 
sphere (oblate or prolate depending on the initial rotation law) 
are enhanced in the magnetized models, and that the MRI is ex- 
pected to develop on the prompt shock propagation time scale. 
ITakiwaki et al.. Q004) observe the formation of a tightly colli- 
mated shock wave along the rotati onal axis for strongly magne- 
tized models. 'Kotake e t al]ll2005l) find an at most 0.5% change 
of the neutrino heating rates ev en in their mos t strong ly mag- 
netized models (\B\ > lO'^GV lArdelian et alJ ( l2005n emplov 
a 2D implicit Lagrangian code, a simplified equation of state, 
and consider energy losses by neutrinos and iron dissociation. 
They add a magnetic field of quadrupole-like symmetry with 
an energy of 10"^ of the core's gravitational binding energy 
to the collapsed, post-bounce differentially rotating, station- 
ary core. The toroidal field component of this seed field first 
grows linearly due to differential rotation, but then starts to am- 
plify exponentially due to the action of the MRI. The resulting 
drastic increase of the magnetic pressure eventual ly causes an 
explosion with an energy of 0.6 10^' erg. Finally. 'Saw a i et alJ 
(>2005i) extend the work of lYamada & Sawai (2004) by consid- 
ering inhomogeneously magnetized cores mainly in the very 
strong field regime {\B\ ~ lO'^G. . . lO'^G), which may pro- 
duce magnetars. They find that poloidal magnetic fields which 



are initially concentrated toward the rotation axis produce more 
energetic explosions and more prolate shocks than cores with 
an initiall y uniform field. A core with an initially quadrupo- 
lar field (kArdelian et al.l2005l) gives rise to a collimated fast jet 
(v ^ c/2), while a core with a pure toroidal fields shows no sign 
of an explosion. 

The paper is organized as follows: we will describe the 
physics included in our models in Sect.|2l and briefly introduce 
our numerical method in Sect.|3 Our results will be discussed 
in Sect.|3 and a summary and conclusions of our work will 
be presented in Sect.|5] In the appendices we provide a brief 
discussion of the relaxing TVD scheme employed in our nu- 
merical code (appendix IaIi. a presentation of the quadrupole 
formula used for the extraction of the GW signal (appendixO. 
and a compilation of some characteristic properties of all mod- 
els (appendixiDt. 

2. Physics of our models 

2.1. Evolution equations 

We evolve the density p, the velocity v, the total energy den- 
sity e^, = e H- eidn + Smag (s, fikin = \pv^ , and e,nag = 
are the internal, kinetic, and magnetic energy density, respec- 
tively), and the magnetic field B of our models using the equa- 
tions of Newtonian ideal magnetohydrodynamics (MHD): 



d,p + v,„(pv"') = 0, 

d, (pv„) + V„, (pv„v"' + - b„b"') = f„, 
d,e^ + V„, ((e* + P,) y'« - W) = q. 



(1) 
(2) 
(3) 



Here, Latin indices run from 1 to 3, and Einstein's sum con- 
vention applies. = Pgas -I- b'^ 12 is the total pressure, which 
is the sum of the gas pressure Pgas and the isotropic magnetic 
pressure P^ag = with b - B/ y/4n. 

Using the ideal MHD equations, we neglect effects due 
to the viscosity and the finite conductivity of the gas. This is 
normally a very good approximation for the stellar interior. 
However, during core collapse interesting hydrodynamic ef- 
fects might arise fr om the inclusion of visc osity, in particular in 
the case of MHD ([Thompson et al.lE004l) . Non-ideal terms in 
the induction equation might lead to reconnection of field lines, 
thus possibly affecting the topology of the field, and might 
prov e important for several kinds of hydromagnetic instabili- 
ties ('Spruit'1999'). 

For a Newtonian self-gravitating fluid the source terms /„ 
and q in the MHD momentum (|5J and energy (0 equations are 
given by 



fn 



-pv VO, 



(4) 
(5) 



where the gravitational potential O obeys the Poisson equation 

A O = AnGp, (6) 

with G being the gravitational constant. The Poisson equa- 
tion is solved in every time step using the solver of 
iMiiller & Steinme'tzl(ll995l) . which is based on the integral form 
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of Poisson's equation and on an expansion of the density dis- 
tribution into spherical harmonics. 

We integrate the MHD equations in spherical coordi- 
nates assuming axisymmetry, i.e. we cannot simulate non- 
axisymmetric instabilities which can occur if the rotation rate 
exceeds a critical value during core co llapse due to angular mo- 
mentum conservation (lTassoullll978i) . Axisymmetry also in- 
hibits the growth of various MHD instabilities (Spruit 1999). 
We further assume equatorial symmetry in order to reduce the 
computational costs of a simulation. 

2.2. Microphysics 

We do neither consider nuclear reactions nor neutrino trans- 
port, and use a simplified equation of state (EOS). Since neu- 
trinos are thought to play a major role in the revival of the 
stalled prompt shock, our approach is limited to the collapse, 
bounce, and shock formation phases when neutrinos are not yet 
dynamically important. This limitation allows us to focus on a 
specific, not yet comprehensively studied part of core collapse, 

namely the influence of MHD eff'ects. 

We have used the approximate, analytic EOS of'janka et al.' 
( I1993I) in our parameter study. It is based on a decomposition 
of the gas pressure P into the sum of a polytropic part Pp and a 
thermal part P±- 



2.3. Initial models 

2.3.1 . Equilibrium models for rotating polytropes 

Concerning rotation and magnetic field, the conditions in the 
stellar core at the onset of collapse are not well constrained. 
Therefore, we have investigated the evolution of a sufficiently 
broad set of simplified stellar models. Except for the magnetic 
field, the initial models are the same as those studied by ZM 
and DFM. 

The initial models are rotating polytro pes in hydrostatic 
equilibrium, constructed with the method of Eriguchi & MiilleJ 
1 19851). They rotate accor ding to the so-called 7-constant law 
I Eriguchi & MiilleJ 19851) 

- ^" (12) 

where the angular velocity Q is given as a function of the cylin- 
drical radial coordinate m. The parameter fio is the maximum 
angular velocity, and the parameter A, having the dimension of 
a length, determines the degree of differential rotation of the 
core. Deviations from rigid body rotation are significant only 
for radii m '» A, where the rotation law approaches that of a 
configuration with constant specific angular momentum j. 



P^Pp + Ah- (7) 
The polytropic part is given by 

Pp^Kp-p^' (8) 
where the adiabatic index 

Fp J^' (9) 
P \ Fa for p > pnuc 

describes cold matter that undergoes a phase transition at nu- 
clear density p^uc = 2 ■ 10'"* cms"'. Below this threshold den- 
sity, the pressure is dominated by a relativistic degenerate elec- 
tron gas with Fp = Fi $ 4/3. At pnuc, the EOS stiffens consid- 
erably, and the adiabatic index jumps to a value Fp = F2 ~ 2.5, 
which mimics the phase transition to incompressible nuclear 
matter. Continuity of the pressure at nuclear density implies 

K2^Ki-pllf\ (10) 

where /ci =4.897 ■ 10''* cgs-units. 

The polytropic pressure changes only due to adiabatic pro- 
cesses, i.e. it cannot describe dissipation of kinetic into thermal 
energy in shocks. This shock heating is treated by the thermal 
part of the EOS which has the form of an ideal gas EOS: 



2.3.2. The hydrodynamic parameter space 

The initial models are calculated using the hybrid EoS in 
the "cold" limit, i.e.Pth = 0, an adiabatic index Fi - 4/3, 
the j -constant rotation law il2\ . and with a central density 
Pc = -10'" cms"'. Instead of f2o we use the ratio of rota- 
tional to gravitational binding energy, /?iot = l^iot/^^giavl (ini- 
tially iigiav ~ 5.5 ■ 10^' erg) to parametrize the models. The 
two-dimensional parameter space {A, Prod is covered by 4 x 5 
initial models according to Tabled 

As Fi =4/3 for our initial models, they are only marginally 
stable against collapse, which can be triggered by either low- 
ering the coefficient ki or the adiabatic index Fi. The for- 
mer approach mimics the energy consumption due to photo- 
disintegration of nuclei, while the reduction of Fi models the 
softening of the EOS due to deleptonisation. We apply the latter 
method to initiate the collapse. 

For realistic equations of state 1.28 ^ Fi 1.325 and 
2.4 ^ F2 ^ 3, respectively. Following ZM and DFM we con- 
structed models with five different values of the sub-nuclear 
adiabatic index Fi, and set F2 = 2.5 (TableQ. The nomencla- 
ture of the models also follows that of ZM and DFM; i.e. model 
A1B3G5 has a rotation parameter A = 5 ■ 10^ cm (Al), a frac- 
tional rotational energy /3 - 0.9% (B3), and a sub-nuclear adi- 
abatic index Fi = 1.28 (G5; see TableQ- 



Ah = (r-i)eth. (11) 

In our simulations, we took y = 1.5 corresponding to a gas 
composed of a mixture of a relativistic and a non-relativistic 
component. 



2.3.3. Magnetic field configuration 

The magnetic fields of our initial models are calculated from 
the vector potential of a circular current loop of radius r,nae in 
the equatorial plane (ii.Iacksoa.1962) . The only non-vanishing 
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Table 1. Initial models and their parametrisation: A and /3rot 
are the rotation law parameter (eauation ll2t and the ratio of ro- 
tational to gravitational energy, respectively. Larger values of A 
correspond to more rigidly rotating cores. Fi is the sub-nuclear 
adiabatic index of our hybrid equation of state (see Sect.|^). 



Model 


A [cm] 


Model 


/3,ot[%] 


Model 


r, 


Al 


5-10^ 


Bl 


~ 0.25 


Gl 


1.325 


A2 


1 • 10« 


B2 


« 0.45 


G2 


1.32 


A3 


5 • 10^ 


B3 


~ 0.9 


G3 


1.31 


A4 


1 • 10^ 


B4 


« 1.8 


G4 


1.30 






B5 


« 4.0 


G5 


1.28 



strengths in iron cores probably do not exceed lO^G, and the 
toroidal field compo nent is expecte d to be much stronger than 
the poloidal one (Heger et alJ2005h . However, as such "weak" 
initial fields do not give rise to important dynamic effects on the 
time scales under consideration here (unless MRl amplification 
would take place; see Sect. 0}, and as we want to investigate 
the principal effects of magnetic fields on the core collapse, we 
consider stronger fields in our parameter study. Weaker initial 
fields may lead to similar effects after a longer amplification 
phase. Our initial fields are purely toroidal, but they rapidly 
(within a fraction of the collapse time scale) develop a strong 
toroidal component. 



component of the magnetic vector potential is A^. It is given by 

A^(r) oc / — d(- cos 0') d(f>' 

sin ff cos 6'5(cos 9'~)6(r — rm„„) 

— —, (13) 



\r-r'\ 

which can be expanded in terms of Legendre-Polynomials f ™, 
yielding 



Z 



.2«+l 



(-l)«(2n- \)\\ri 
2"(n+ l)!r2«+2 

= min(r, r^„) and r 



n,+i(cos0). 



(14) 



Here, r< = min(r, r^ag) and r> - max(r, r^ag) and (2n - 1)!! 
= 1 ■ 3 ■ 5 ■ ... ■ (2n + 1). The constant of proportionality in 
these expressions is fixed by demanding that the magnetic field 
strength in the core is equal to a given value. 

The magnetic field strength in the core's center is normal- 
ized by Bo = V47r/7o- For very small radii (r <k r^ag) the field 
resembles a uniform field parallel to the rotation axis, whereas 
at large radii the field lines bend towards the equatorial plane. 
The field strength is largest in the interior of the current loop 
(Fig.Q. 

Table|2provides an overview of the parametrisation of the 
magnetic field. The models are denoted as follows: Model 
A3B3G3-D3M12 is the hydrodynamic initial model A3B3G3 
(Sec . 12.3.21 endowed with a magnetic field generated by a cur- 
rent loop located at a radius of r^ag = 400 km and a maximum 
field strength of bo = 10'^ G (see Tabled}- The initial magnetic 
energy E^^g for models AaBbGg-DdM12 is given in Table|3] 

Most simulations were performed using models AaBbGg- 
D3Mm (m = 10,11,12,13). Their magnetic field config- 
uration is shown in Fig.^ Since the radius of the current 
loop (rniag = 400 km) that generates this field configuration 
is small compared to the radius of the stellar core (rcore ~ 
1500... 2 100 km), the magnetic energy is highly concen- 
trated in the center of the stellar core. This is different for 
models AaBbGg-DOMm, which posses a homogeneous field 
directed along the rotational axis. The magnetic energy of 
these "uniform-field" models is much larger than that of the 
corresponding "current-loop" models AaBbGg-DdMm (d = 
1,2,3,4) due to the contributions of the outer layers of the core. 

The magnetic field strengths of our initial models are - as in 
most studies of magneto-rotational collapse - much higher than 
those estimated to exist in realistic stellar cores. Magnetic field 



Table 2. Parametrisation of the initial magnetic fields for 
the models of series AaBbGg-DdMm by the radius of the 
field generating current loop centered at r^ag (parametrized 
by c/ = 1,2,3,4,0) and the field strength in the core's center 
Bq - V47r/7o (parametrized by m - 10, 1 1, 12, 13). For models 
AaBbGg-DOMm the field generating current loop is located at 
infinity yielding a uniform magnetic field throughout the entire 
core. 



Model 


'■mag [km] 


Model 


bolG] 


Dl 


100 


MIO 


10'" 


D2 


200 


Mil 


10" 


D3 


400 


M12 


10'2 


D4 


800 


M13 


10'3 


DO 


oo 







Table 3. Initial magnetic energy /Smag and typical values of the 
ratio of magnetic to gravitational energy jSmag (the exact values 
depend also on the hydrodynamic initial model and its grav- 
itational energy) for the models of series AaBbGg-DdM12. 
The magnetic energy of models with bo + 10 '^G can be 
obtained by a simple scaling relation, e.g. £'mag^''^ oimw _ 



Model 


£'mag [erg] 


log/3mag 


DOM 12 


1 ■ 10^^^ 


-2.7 


D4M12 


6 ■ lO^** 


-3.0 


D3M12 


7 • 10*^ 


-3.9 


D2M12 


7 • 10^"^ 


-4.9 


D1M12 


1 • lO'"^ 


-5.8 



2.4. Gravitational-wave emission 

During collapse, bounce, and explosion the rapid infall of mat- 
ter and in particular its more or less abrupt slowdown give rise 
to strong variations of the matter-density quadrupole moment 
of any aspheric core. This causes the emission of gravitational 
radiation. 

We calculate the gravitational wave amplitude of the core 
using the quadrupole formula in spherical coordinates , and ap- 
plying the extension of the formulation of iMonchmever et alJ 
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28.32 



Fig. 1. The initial magnetic field configuration of the models 
of series AaBbGg-D3M13. Besides the field lines the distribu- 
tion of the magnetic pressure (color coded) is displayed. P„iag 
is largest near the center of the field generating current loop lo- 
cated at a radius of r^l^ = 400 km in the equatorial plane, and 
drops rapidly with increasing radius for r > r™g. 



|l99l[ MSMK. hereafter) to the MHD case due to lKotake et alJ 
( 1200 4 m. Our treatment includes the hydrodynamic, gravita- 
tional, and magnetic forces acting on the fluid. We calculate 
the quadrupole amplitude Ajq according to the formula (see 
AppendixO: 



4E2 _ 
^20 - 



G 32712 

VTf 



Jo Jo 3 



frr(3z^-l) + feg(2-3z')-f^^ 



where the components of fij are given by 

fij = pViVj - bibj. 



(15) 



(16) 



In the following, we will refer to the various parts of the to- 



tal amplitude as follows: A^q.^, ^, 



aE2 



and AjQ.j, denote 



the contributions of the terms involving v,Vj, bibj, and 5,0 in 
Eq.[T5l respectively. Furthermore AH^^^, A^.^^^, and A^.^^^^ 
are the sums over all components of A 
respectively. 



^20 Ai, ' ^"'^ ^20;G, ' 



E2 

20;\vi', 



The radiative quadrupole moment Mjq (Ea. lC.3> is a mea- 
sure of the asphericity of the core's density distribution. It is 
positive for a very prolate core, and negative in the limit of very 
oblate cores. Its first time derivative A^jo (Eq-1^3} measures 
the asphericity of the mass-flux and the momentum distribu- 
tion of the core, and its second time derivative, the quadrupole 



amplitude A^q , is a measure of the asphericity of the forces act- 
ing on the fluid. As a rule of thumb, a prolate mass-flux or a 
prolate momentum distribution (e.g. a bipolar jet-like outflow 
along the rotational axis) gives rise to a positive value of A^^o^. 
Forces that act on the core in a way to make it more oblate, 
such as the centrifugal force that has its manifestation in the 
part of the amplitude, will give rise to a negative contri- 



^E2 
20;i>iv 



E2 

2O;ivi'0 



term 



bution to the total amplitude (negative sign of the A 
in Eq.[l5}. 

The different signs of the hydrodynamic and the magnetic 
contributions to the amplitude (Ea.ll6> resulting from the dif- 
ferent signs of the hydrodynamic (Reynolds) and the magnetic 
(Maxwell) stresses in the MHD flux terms, will - for suited 
topologies of field and flow - lead to a more or less prominent 
phase shift between the hydrodynamic and the magnetic am- 
plitude. If the gravitational wave amplitude A^,^ is in phase 
with the hydrodynamic amplitude A^y^^^ (which holds well for 
many models, in particular for those with relatively long oscil- 
lation periods; see Sect.|^, the magnetic amplitude A^y^,^^ may 
be phase shifted with respect to A^,?.^,,^ + A^^.r- ■ Such a phase 
shift was observed bv .Yamada & Sawai. C2004.) . 



3. Numerical method 

The MHD equations are integrated using a newly developed 
Eulerian, finite volume code based on the algorithm devised by 
iPen et al. ( 2003). Th is code employs the relaxing TVD method 



of Jin & Xig ( 1 19951) for the solution of the ad vection equations 
and the constraint-transport formulation of lEvans & HawlevI 
fl988.) to deal with the divergence constraint of the magnetic 
field. 

We have rewritten the original code of Pen et al. in order 
to adjust for the simulations of stellar core collapse. This in- 
cluded the transformation of the equations from Cartesian to 
spherical coordinates, the calculation of the gravitational po- 
tential, the implementation of the gravitational source terms in 
the momentum and energy equations, and the implementation 
of an approximate equation of state for iron core matter. The 
integration of the fluid equations and of the induction equation 
is based on a second-order (piecewise-linear) relaxing TVD 
method. For a short summary of this method see AppendixlAl 
For the time evolutio n we use an oper ator-split approach based 
on a method of lines (lLeVeauell992h . 

The simulations were performed on a grid of 380 loga- 
rithmically spaced radial zones up to 7?* km, where 7?* ~ 
1700... 2000 km is the radius of the initial stellar model. The 
central resolution was (A,)c ~ 300 m. The angular grid con- 
sisted of 60 equidistant zones in the domain < 6* < |. This 
grid resolution has been chosen after obtaining converged re- 
sults when running several models at different resolutions. The 
numerical convergence of our simulations is demonstrated in 
the Appendix (Sect.|B}. 

The stellar models used by us (polytropes) are quite com- 
pact configurations of matter characterized by a sharp transition 
from a high density interior to a low density surface layer. Note 
that this feature is also seen in sophisticated stellar evolution 
calculations, which predict a steep density gradient at the outer 
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edge of the iron core. In rotating models the transition layer is 
aspherical and must completely be contained within the spher- 
ical boundary of the numerical grid, i.e. its numerical treatment 
requires special care. Grid zones outside the core are filled with 
an "atmosphere" fluid of some prescribed density Patmo 31 rest. 
During collapse the infall of matter creates a region near the 
edge of the core where the density might become so low that 
numerical problems arise. To overcome these difficulties, we 
follow the approach of DFMl and set the hydrodynamic vari- 
ables p, v,e equal to some prescribed values patmo, 0, Catmo in 
all zones where the fluid density falls short of a given thresh- 
old pcut- In this way the atmosphere can adjust to the (non- 
spherical) varying shape of the star. We used pcut = 10^ cm s"', 
and set Patmo = 10^ cms From the atmospheric density and 
velocity one can calculate the energy density assuming zero 
thermal energy. 

The evolution of the magnetic field is turned off for zones 
marked as atmosphere, i.e. b remains constant in the atmo- 
sphere consistent with the assumed zero velocity of the atmo- 
sphere gas. 

4. Results 

The results of our simulations show that magneto-rotational 
collapse can be categorized in essence into two limiting cases 
depending on the strength of the initial magnetic field. If the 
initial magnetic field is weak, its influence on the dynamics 
and the gravitational wave emission is negligible during the 
time scales of our simulations (see Sect. l4.H . The results of the 
hydrodynamic simulations of ZM and DFM apply in this case 
without any modification. On the other hand, and not unex- 
pectedly, initially strongly magnetized cores evolve quite dif- 
ferently, as will be discussed in Sect. l4.2l Note that since the 
MRI acts independently of the strength of the initial magnetic 
field, the distinction between dynamically negligible and dy- 
namically important fields may be an artificial one resulting 
from our inability to simulate the MRI for magnetic fields be- 
low a certain threshold (see Sect. l4.3> . Hence, also initially 
weak magnetic fields may cause similar dynamical effects as 
strong ones. The dependence of our results on the initial mag- 
netic field configuration will be discussed in Sect. l4.4l and fur- 
ther information about the temporal evolution of all models is 
provided in appendixiDl 

4.1. Weak initial fields 

4.1 .1 . GW signal and dynamics 

The gravitational wave signals of initiafly weakly magnetized 
cores do not differ from those obtained by MM, ZM and DFM 
for the corresponding non-magnetized initial models, because 
the magnetic fields never become dynamically important dur- 
ing the simulations. Consequently, both the magnetic force 
contribution to the GW amplitude (see Sect. l2.4l i and the back- 
reaction of the magnetic field on the flow are negligible. This 
"weak-field" behavior holds for most models with initial fields 
of ^ 10" G. However, in some models bouncing at relative 
low central density due to centrifugally forces (e.g. models 



A4B5G5-D3M12 and A2B4G1-D3M12) even a field of 10'^ G 
does not affect the evolution of the core significantly. We note 
that these numbers may change when the field amplification by 
the MRI, which works independently of the seed field strength, 
is properly simulated (see Sect. 14. 3> . 

Three types of GW signals are distinguished by MM 
and ZM which they call type-I, type-II, and type-Ill, respec- 
tively. The evolution of three representative weak-field mod- 
els that emit these different signal types is illustrated in Fig.|2l 
Although the GW signals are very similar to those of the coiTe- 
sponding non-magnetized cores of ZM, we will discuss them 
and the underlying dynamics here in some detail for reason of 
comparison with the GW signals of strongly magnetized cores 
(see Sect.lO. 

In the case of a type-I (or standard-type) GW signal (model 
A1B3G3-D3M10; Fig.|2 left panels), the core bounces due to 
the stiffening of the nuclear EOS at supra-nuclear densities. 
Pressure waves crossing the inner core stop the infall and lead 
to the formation of an outward moving shock wave. Hence, the 
typical time scale is roughly given by the sound crossing time 
of the inner core at bounce (~ 1 ms). After bounce the core 
exhibits damped oscillations on roughly the same time scale. 
During collapse, the GW amplitude is increasingly positive due 
to the gravity contribution Ajq.^^.^^, > which dominates the 
hydrodynamic one AjQ.jjy^ < 0. At bounce the GW amplitude 
decreases strongly assuming large negative values, because the 
modulus of the centrifugal contribution Ajg.j,^,,^ becomes very 
large when the rotational energy approaches its maximum. The 
oscillations of the core produce the oscillations of AS? and 

^ 2U;grav 

hence of the total GW amplitude on the same time scale. As 
the shock wave is almost spherical (except for initially very 
rapidly rotating cores), the post-bounce GW amplitude is pre- 
dominantly produced by the central core with only a modest 
contribution of the outer layers. The hydrodynamic contribu- 
tion A5?. , is dominated by A?,? 

2();hyd zlnv^v^ 

A type-II signal (model A2B4G1-D3M10; Fig.|3 middle 
panels) is emitted, if, for an initially sufficiently fast rotation 
and a sufficiently large degree of differential rotation, the core 
suffers a bounce due to centrifugal forces before nuclear den- 
sity is reached. After bounce the core exhibits several large am- 
plitude oscillations, which are considerably less damped than 
in the case of a pressure dominated bounce. The oscillations 
are reflected in the GW signal: each time the density reaches a 
maximum the GW amplitude becomes strongly negative, while 
being positive otherwise. The negative GW amplitude again re- 
sults from the enhanced rotational energy during this collapse 
phase when the inner core also becomes strongly oblate. The 
contribution of the centrifugal-force amplitude A^^y^,^^,^ to the 

hydrodynamic one Ajq ^^^^j is negligible. Since much of the dy- 
namics happens at relatively large radii, the GW signal is pro- 
duced by the whole core. The contribution of the central core 
is typically negative because of its rotationally flattened, oblate 
shape, whereas the outer layers contribute to the overall am- 
plitude with a positive signal due to the prolate shape of the 
outward propagating shock wave. 

The signals of type III (model A3B3G5-D3M10; Fig.H 
right panels) that are emitted by cores collapsing very rapidly 




Fig. 2. The evolution of the maximum density (top panels) and the GW amplitude (bottom panels) of three weak-field models: 
the standai-d type-I model A1B3G3-D3M10 (left panels), the type-II model A2B4G1-D3M10 (middle panels), and the type-Ill 
model A3B3G5-D3M10 (right panels). The bottom panels show the total GW amplitude Ajg (solid lines) together with the partial 
amplitudes AjQ.fjyj (dashed lines), A^q.^^^^ (dash-dotted lines; here almost zero), and A^q.^^.^^ (dash-dot-dot-dot line). 



due to a very soft sub-nuclear EOS are characterized by the 
lack of a strongly negative GW amplitude at bounce. This 
results from the only modest rise of the rotational energy at 
bounce implying a relatively small centrifugal contribution to 
the GW signal A^q.^,^^,^, and from the relatively large positive 
contribution of the still collapsing layers outside to the shock 
wave. Past bounce the signal exhibits rapid (on the hydrody- 
namic time scale fhyd ~ l,ms) oscillations like the maximum 
density. The contributions from the outer layers of the core are 
also responsible that the signal then remains positive for an ex- 
tended period of time (i.e. several oscillation periods) during 
which the amplitude gradually decreases. 



4.1.2. Field amplification 

The evolution of the magnetic energy is illustrated for four se- 
lected models in Fig.|3]showing the ratio of the field energy and 
the gravitational energy for the total iPmug), toroidal (J3^), and 
poloidal (fisnan - P<ti) magnetic field, respectively. Both the total 
and the toroidal magnetic energies rise sharply at bounce by a 
factor of 10^ to 10^. The magnetic energy is mostly stored in 
form of the toroidal magnetic field B^, which is created from 
the poloidal component by the action of differential rotation. 



The amplification process extends beyond the end of our simu- 
lations. 

As our study is restricted to axisymmetric simulations, the 
main field amplification mechanism is the conversion of rota- 
tional into magnetic energy via the Q-dynamo. Axisymmetry 
suppresses most of the instabilities of the (toroidal) field that 
are necessary to close the dynamo loop of a full a-Q. dynamo, 
where a poloidal field is converted into a toroidal one and am- 
plified by differential rotation, and then converted back into a 
poloidal fie ld by some instability of the toroidal field dSpruitl 
Il999ll2002l) . 

Before we discuss the properties of some models in greater 
detail, we summarize a few general trends: 

- The larger the initial rotational energy (for a given degree 
of differential rotation and a given EOS), the larger is the 
rate of field amplification. 

- The higher the degree of differential rotation (for a given 
Plot), the larger is the rate of field amplification. 

- Among a series of models with the same initial configu- 
ration, the magnetic field is amplified more efficiently for 
models whose EOS has a larger sub-nuclear adiabatic in- 
dex Fi and do not suffer a centrifugal bounce. The time 
scale for the amplification of the magnetic field is set by 
the rotation period of the core. 
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Fig. 3. The evolution of the magnetic energy parameters /S^ag 
(total energy, solid line), /S^ (toroidal field energy, dashed 
line), and jSmag - (poloidal field energy, dotted line) of 
models A1B3G3-D3M10 (upper left panel), A3B3G5-D3M10 
(upper right panel), A2B4G1-D3M10 (lower left panel), and 
A4B5G5-D3M10 (lower right panel), respectively. Note that in 
the upper left panel /3m-dg ~ yS^ and, thus, both lines are almost 
indistinguishable. 

- The collapse of models with a relatively stiff sub-nuclear 
EOS (AaBbGl-DdMm) leads to a very extended post- 
bounce core having a rotational frequency much less than 
that of a coiTesponding type-I model even though the ro- 
tational energy and /3roi may be larger in the former model. 
Thus, the field is amplified more slowly in models suffering 
a sub-nuclear bounce due to centrifugal forces, and their 
cores experience several phases of expansion during which 
they slow down. This leads to a significantly less efficient 
amplification or even a weakening of the magnetic field. 
We thus find that type-I models are most efficient in terms 
of magnetic field amplification. 

- Due to the efficient field amplification most of our weak- 
field models reach magnetic field strengths of the order of 
lO'^G to lO^^G within a few tens of milUseconds after core 
bounce. This is in the range of field strengths observed in 
magnetars. From the definition of jSmag follows, that 

where bis a typical value of the field within the core, and R 
and M are its mass and radius, respectively. Thus, a typical 
post-bounce field corresponds to a magnetic energy param- 
eter of 

^™.~0.05%(^[(^f(^)\ (18) 

The maximum energy that can go into the magnetic field is 
limited by the energy that is contained in the differential rota- 



tion of the core. None of our weak-field cores is evolved long 
enough for the field to reach saturation strength. However, ex- 
amining the field amplification process for models with strong 
(^ lO'^G) initial fields, which reach dynamically important 
values not too early, we conclude that efficient field amplifica- 
tion ceases when the magnetic energy of the collapsed core ap- 
proaches a significant fraction of its rotational energy, i.e. some 
ten percent. As all of our cores initially have similar rotational 
energies (/6,ot about a few percent), we expect them to reach 
(surface) field strengths of the same order if there is sufficient 
time for amplification. From our models with initially stronger 
fields we can estimate the saturation field strengths to be of the 
order of b ^ 10'^ G. Using the estimate for /Jmag given above 
( I18> . these fields correspond to jSmac; ~ (0.05. ..5) %. Such val- 
ues we have found in our simulations. The magnetic field in 
most post-bounce cores is strongly concentrated in the inner 
core, reflecting the density structure. Thus, for usual standard- 
type cores, the field is strongest a few kilometres from the cen- 
ter, and drops rapidly towards larger radii. When a strong shear 
layer is present at the surface of the inner core, the field may 
also be strongest there. This is the case for many of the weak- 
field cores, but does not necessarily hold for the strong-field 
ones. In the latter case, the simulations tend to give values of 
ySmag in excess of the previous estimate. 

The field amplification via the Q-dynamo is most efficient, 
if the poloidal field component is large compared to the toroidal 
one. A characteristic growth time for the generation of b^ from 
a radial component bj;j is given by (Meier et al..l976i) 



bj^Tudj^Q. ' 

In our models the poloidal component does not grow after 
bounce (apart from the effects of compression and expansion 
during the large-scale oscillations of type-II models) unless 
there is an instability (Fig. 13 Sect. l4.3> . 

The rate of field amplification is similar for the fast collaps- 
ing cores A1B3G5-D3M10, A2B4G5-D3M10 and A3B3G5- 
D3M10, which have a soft sub-nuclear EOS and which do not 
suffer a centrifugal bounce. A slightly more efficient amplifica- 
tion is observed for the initially more differentially and faster 
rotating cores A2B4G5-D3M10 and A3B3G5-D3M10 than for 
the initially rigidly rotating model A1B3G5-D3M10. An even 
faster amplification is observed for the slower collapsing model 
A1B3G3-D3M10, where 10 ms after bounce ySmag = 8 ■ 10"^ 
compared to y6n,ag = 8-10"''formodel A3B3G5-D3M10(Fig.|3l 
upper two panels). This reflects the different infall profiles of 
the cores whereby the regions of strongest magnetic field (ini- 
tially around and interior to the field generating cuiTent loop 
at r^iag - 400 km) are swept along to different positions in the 
core. In case of the model A1B3G3-D3M10 this region is much 
closer (r » 40 km) to the center at the time of bounce than for 
model A3B3G5-D3M10, where it is located at r a; 120 km. 
Therefore, the fraction of the magnetic field that is available in 
the regions where the Q-dynamo acts most efficiently is larger 
for the former model, i.e. the field amphfication can proceed 
more efficiently. 

During collapse and to an even larger extent during and 
after core bounce the magnetic field configuration of model 
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Fig. 4. Magnetic pressure (color shaded), velocity field (ar- 
rows), and magnetic field lines of model A1B3G3-D3M10 at 
t = 51.44 ms. 
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Fig. 5. Same as Fig.Hl but for model A4B5G5-D3M10 at f 
36.13 ms. 



A1B3G3-D3M10 is considerably distorted. The field lines ini- 
tially surrounding the off-center current loop (Fig.[3 are pulled 
towards the center. Apart from an overall compression the field 
geometry remains basically the same as in the initial model. 
A complex pattern of "filamentary" regions of high and low 
magnetic fields develops as the field lines get entangled by the 
fluid flow around and after bounce. No simple classification of 



the magnetic field as a dipole, quadrupole, etc., is possible. In 
the subsequent evolution the tangled fields follow the outward 
propagating shock front which is almost spherically symmetric 
in this initially rigidly rotating model (Fig.0}. 

Multiple-bounce models that experience several phases of 
collapse and expansion exhibit a similar evolution of the mag- 
netic field. During a contraction phase angular momentum con- 
servation yields a more efficient amplification of the magnetic 
field energy, whereas ySmag grows much slower in an expan- 
sion phase. This is the case for model A2B4G1-D3M10 (Fig.E 
lower left panel). The magnetic energy parameter /Jmag rises by 
a factor ^ 3 during 10 ms around core bounce {% x 100 ms), 
but then requires 40 ms to amplify the field by another factor 
of a: 1 .6. The amplification rate increases again strongly during 
the subsequent contraction phase. For model A4B5G5-D3M10, 
the amplification factor during core bounce is about 20, but 
later the field energy gets greatly reduced (/Sniag drops by a fac- 
tor a! 3.5) during the extremely violent expansion that follows 
the bounce in this model (Fig.|3] lower right panel; the evolu- 
tion of the maximum density is shown in Fig.ll9>. Afterwards, 
differential rotation again starts to increase the magnetic field, 
but on a longer time scale. 

Model A4B5G5-D3M10 possesses a torus shaped initial 
density distribution, and maintains it throughout the entire 
evolution. The field structure of this model tends to become 
very complex during the evolution. In the final model, at t = 
36.13 ms, the field exhibits both sheet-like regions of strong 
fields (a 100 km off the axis), and also toroidal (at r a 120 km, 
z ~ 60 km) and cylindrical (e.g. a 50 km off the axis) weak- 
field regions (Fig.|5}- 

4.2. The Strong field case 

For initially strong magnetic fields the collapse shows signifi- 
cant deviations from the purely hydrodynamic case. The most 
striking new feature is the braking of the rotation of the core 
by magnetic stresses. The initially strongest fields are ampli- 
fied by differential rotation during collapse to such a level that 
they can cause considerable braking, whereas initial fields of 
the order of 10'^ G or less require the action of an MRl-like 
instability to reach a level sufficient for significant angular mo- 
mentum transport on the collapse time scale. In the following 
we will first discuss the former class of models divided into 
models bouncing due to pressure forces ('Sect. l4.2"Tl or cen- 
trifugal forces (Sect. l4.2'!2t . and later address the issue of the 
MRl in Sect.lO 

4.2.1 . Models bouncing due to pressure forces 

As in the case of weak-field models, the poloidal field energy 
of a single bounce models such as A1B3G3-D3M12 is approx- 
imately constant during the early post-bounce evolution, but 
it is strongly amplified afterwards by MRI-like modes in the 
core which efficiently remove angular momentum from the in- 
ner core and decrease its rotational energy. The evolution of 
the rotational energy hence shows a clear deviation from that 
of a weak-field model (Fig.|6}, which however should not be 
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Fig. 6. The dynamics and the GW signal of model A1B3G3-D3M12. The upper left and middle panels provide a comparison 
of the evolution of the rotational energy parameter jSrot and the maximum density of this model (solid lines) with those of model 
A1B3G3-D3M10 (dashed lines). The upper right panel displays the evolution of jSmag (solid line), yS^ (dashed line), and /Jmag - 
(dotted line), respectively. The GW signal is shown in the lower panels: total amplitude (solid line, left; the dashed line gives the 
total ampHtude of model A1B3G3-D3M10), -A\l^^^^ (solid Hne, middle), A^l^^.^^ (dashed Une, middle), and A^^^^^ (right). The 
corresponding variables of model A1B3G3-D3M10 are shown in Fig.|2 



the case, as the MRl does not depend on the seed field. We will 
explain in Sect. l4.3l whv this different behavior arises. 

The loss of rotational support also affects the density struc- 
ture. The core of model A1B3G3-D3M12 is slightly less ro- 
tationally flattened than the core of a non-magnetized model 
such as A1B3G3-D3M10. The shock that forms at bounce, 
when the models hardly differ, is unaffected by the different 
evolution of the central region, whereas matter in the central 
region experiences an additional acceleration due to the pres- 
ence of the strong magnetic field. The decrease of the rotation 
rate, and also of the degree of differential rotation in the strong- 
field case limits the amplification of the toroidal field compo- 
nent by wrapping of field lines (Fig.|6l upper right panel), and 
the total magnetic energy has a significant poloidal contribu- 
tion. At f a: 70 ms, the ratio of toroidal and total magnetic 
energy has dropped to a value of about 80 %. The poloidal 
component grows a little during the post-bounce evolution by 
the action of meridional motions in the core. When the mag- 
netic field becomes dynamically important, magnetic energy 
is used to accelerate the fluid, which puts an end to the phase 
of efficient field amplification. During the late stages of evo- 
lution, the total field energy remains approximately constant. 
Its final value corresponds to ySmag ~ 1 %. The field structure 



of model A1B3G3-D3M12 differs from that of the weak-field 
core A1B3G3-D3M10 in the central parts of the core at late 
times due to the MRl modes (that could not be resolved for the 
weaker field; see Sect. l4.3> . but farther away from the center 
the fields have a similar topology. 

Around the time of bounce the GW signal of model 
A1B3G3-D3M12 is very similar to the one of the less magne- 
tized core A1B3G3-D3M10, and also the contributions of the 
individual parts of the total signal are comparable (Fig.|6j lower 
panels). In the post-bounce evolution, both the hydrodynamic 
and the gravitational part of the amplitude start to decrease in 
magnitude, as the rotationally induced asphericity of the core 
decreases along with the extraction of rotational energy. The 
magnetic amplitude increases strongly as magnetic forces act 
on the core decreasing its rotation. The oscillations of the core 
are still imprinted on the post-bounce GW amplitude, but their 
impact on the signal decreases with time. In the long run, the 
amplitude assumes positive values and varies relatively little. 
Dynamically, this phase is characterized by the emergence of 
a weak outflow along the polar axis far behind the shock wave 
created at core bounce. This outflow is predominantly driven 
by magnetic forces (Fig.Q. The positive long-term GW am- 
plitude is mainly due to the Ajq.^, ^ term which is large inside 
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Fig. 7. Snapshot of the inner regions of model A1B3G3- 
D3M12 at f = 76.94 ms showing the magnetic pressure (color 
shaded), the bipolar flow field (vectors), and the magnetic field 
lines, respectively. 



the outflow, whereas the inner parts of the model and the al- 
most spherical shock wave contribute little to the total signal. 
Despite the differences in the post-bounce GW amplitude con- 
tributions of the two models A1B3G3-D3M10 and A1B3G3- 
D3M12 their total signals are quite similar, i.e. distinguishing 
the two models observationally would be very difficult. 

Similar features as those observed in the case of model 
A1B3G3-D3M12 (Fig.|6j are present in the evolution of 
the stronger magnetized model A1B3G3-D3M13 (Fig.|8ll, but 
they are much more prominent. Unlike all previously dis- 
cussed models, this core shows significant deviations from the 
purely hydrodynamic case already during collapse. For model 
A1B3G3-D3M13 the core bounce is delayed by ^ 1 ms com- 
pared to the weaker magnetized models A1B3G3-D3M10...12. 
The rotational energy reaches a maximum value of jS™" s; 7 % 
near bounce, which is lower than in the less magnetized models 
08™" ~ 8 %). Within the next ^ 6 ms the core looses ^ 80 % of 
the rotational energy it had acquired at bounce. Subsequently, 
the rate of energy loss decreases, and the core contracts (Fig.|8j. 
Its maximum density rises from pmax = 2.7 • lO^'* cm s"' in the 
first post-bounce density minimum to pmax = 3.2 ■ 10^"^ cm s"' 
at f = 58 ms before it starts to decreases by a small amount 
again. The process of magnetic braking prevents a further am- 
plification of the field, and it prevents the transformation of 
poloidal into toroidal magnetic field. The ratio jS^/ySmag does 
not approach unity, as in case of the models with weaker initial 
fields, but it remains below ^60% decreasing slowly later. At 
t X 71ms, when we stopped the computations, the value of jSmag 
had dropped to jSmag ~ 2.05 %, and we expect it to eventually 
decrease to a value of ~ 1 %. 



In model A1B3G3-D3M12 the immediate post-bounce 
core is very similar to the corresponding non-magnetic one as 
the magnetic field is not sufficiently strong to change the dy- 
namics unless an MRl-like instability sets in. This is different 
from the more strongly magnetized model A1B3G3-D3M13, 
where without the MRl the combined amplification due to con- 
traction and differential rotation is sufficient to cause the brak- 
ing effect. This is generally true for all of our models: very 
strong initial fields (> 10'^ G) do not need amplification by 
an instability to brake the core, whereas weaker ones heavily 
depend on it. 

After bounce, the shock wave develops a non-spherical, 
bulb-like shape , and the twisted field lines give rise to a highly 
magnetized post-shock fluid that pushes the shock near the 
equatorial plane around 6 70° (Fig.|5). Near the axis, the 
immediate post-shock matter is only weakly magnetized, the 
ratio of gas pressure and magnetic pressure being largest at a 
significantly smaller radius (rmax.mag ~ 27 km) than the shock 
position (rshock ~ 40 km). The shock maintains its bulb-like 
shape also during the secondary contraction phase that goes 
along with the extraction of rotational energy, and fast outflow 
both at high and low latitudes. At later times {t - 60.40 ms, 
Fig.llO>. the fluid pattern has changed. The formerly impor- 
tant outflow near the equator stalls, and the shock surface be- 
comes more elongated. The near-axis outflow velocity has fur- 
ther increased, its maximum value being v™^" - 0.19c in the 
strongly magnetized region well behind the shock (rmaxmag = 
270 km < Tshock = 330 km). In the outflow along the rota- 
tional axis, a cylindrical shell of very low pressure gas can be 
identified (at 6 « 10° in Fig.llO>. The structure of the mag- 
netic field is dominated by two extended highly magnetized 
regions, one located near the equator at (r ^ 200 km, 6 ^ 80°), 
and a cylindrical one oriented along the rotational axis having 
a width of ^ 100km (left panel of Fig.llO>. At an even later 
epoch, f - 68.51 ms (Fig. !! 1> . only little has changed near the 
equatorial plane, however at the axis the highly magnetic gas 
is about to catch up with the shock wave, which is now signif- 
icantly prolate with an axis ratio of about 3 : 2. The outflow 
has accelerated further (v^ = 8 ■ 10^ cms"' ~ c/4), the maxi- 
mum velocity corresponding to the highly magnetized regions 
behind the shock at r = 870 km. Hence, we see the formation 
of a jet-like outflow from the core. 

In its very interior the core is rotating slower than its non- 
magnetic variant at late times. Furthermore, it has developed 
a roughly toroidal region of relatively slow and highly rigid 
counter-rotation. This region grows from a layer at r » 10 km, 
i.e. near the edge of the inner core, which experiences suffi- 
ciently strong magnetic stresses to reverse its direction of ro- 
tation. Interior to the retrograde rotating region the angular ve- 
locity decreases with time. 

Assuming that the total energy density of a fluid element, 
Cfoi - s + jpv^ + jb^ - p ■ O, is converted entirely into ki- 
netic energy ^pv^, terminal outflow velocities Voo of up to 
1.85 ■ 10"' ems ' are observed for model A1B3G3-D3M13 at 
late times, particularly in the bipolar outflow along the rotation 
axis (Fig.^3. The large velocities stretch our non-relativistic 
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Fig. 8. The dynamical evolution and the GW signal of model A1B3G3-D3M13. The upper panels display the temporal evolution 
of the rotational energy parameter y6,ot (left), the maximum density p,nax (middle), and the magnetic energy parameters (right; 
jSmag (solid line), fi^ (dashed line), and jSmag - y60 (dotted line)), respectively. The GW signal of the model is shown in the lower 
panels: total amplitude (solid line, left; the dashed line gives the contribution of the layers with r < 59.4 km), -Ajo.fj^j (solid line, 
middle), J^.^^^^ (dashed hne, middle), and Ajq.^^^ (right). 



MHD approach, and imply that a realistic simulation of the late 
evolution of the outflow will require relativistic MHD. 

The GW ampHtude of model A1B3G3-D3M13 (Fig.|8] 
lower panels) is enhanced by about 10 % at core bounce. After 
bounce, ringing on the dynamical time scale of the core the av- 
erage amplitude shifts to positive values. The amplitudes and 
frequencies of the oscillations are hardly affected by the over- 
all shift of the average amplitude to positive values implying a 
different origin of the oscillations and the overall shift. A few 
milliseconds after bounce grows to ~ 400 cm, and the rapid 
variations of the signal cease. The hydrodynamic and gravita- 
tional contributions to the total amplitude are greatly reduced 
in magnitude. The inner core has lost a large amount of its rota- 
tional energy and is of almost spherical shape. Its contribution 
to the total signal amplitude is hence relatively small for epochs 
well after bounce. 

While the hydrodynamic amplitude produced by the inner- 
most layers of the core is very small, the gravitational and mag- 
netic contributions reach quite large values of several 10^ cm. 
However, as both contributions are of opposite sign they cancel 
each other almost completely (Fig.|8|i, i.e. only a small net am- 
plitude results. The relative smallness of the hydrodynamic am- 
phtude of the central core indicates that - once most of the rota- 



tional energy is extracted - magnetic forces may become more 
important for the core's structure than the centrifugal ones. 

The very non-spherical shape of the shock wave, and in 
particular the appearance of the jet-like outflow that c allies a 
large (radial) kinetic energy, causes a long-lasting positive GW 
amplitude the major contribution being the Ajg.^, amplitude of 
the outflow (Fig. ll2l i. Growing larger with time this contribu- 
tion eventually makes the hydrodynamic amplitude to become 
positive at f = 6 1 .5 ms. Thus, a bipolar motion can be identified 
in the GW signature through the appearance of a long-lasting 
positive signal. 

The behavior just described holds for models A1B3G3/5- 
D3Mm, A2B4G5-D3Mm, and A3B3G5-D3Mm, too. They dif- 
fer, however, concerning the time scales and the vigorousness 
of the phenomena. Most dramatic is the evolution of model 
A3B3G5-D3M13 (Fig.|8}, where the loss of rotational energy 
allows the core to contract to densities (p^ax - 4.4- 10''* cm s"') 
that exceed the bounce density pb - 3.46 ■ lO''* cms ' by up 
to 30 %. The shock wave is already strongly prolate when it 
forms, and the magnetic field is highly concentrated both to- 
wards the equator and - like in model A1B3G3-D3M13 - in a 
cylindrical region oriented along the rotation axis. The rapidly 
moving, highly magnetized outflow deforms the shock wave 
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Fig. 9. Snapshot of model A1B3G3-D3M13 att - 50.19 ms, when the maximum density reaches the first post-bounce minimum. 
In the left panel the flow field (vectors) and the magnetic field (field lines) are displayed together with the magnetic pressure 
(color-coded). The right panel shows the color-coded velocities of in-flowing (i.e. Vr > 0) and out-flowing matter (i.e. v, < 0), 
respectively. 




Fig. 10. Same as Fig.|51 but at t - 60.40 ms. Note that the length scales are changed compared to Fig.|51 
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Fig. 11. Same as Fig.|9] but at f = 68.51 ms. Note that the length scales of the plot are changed compared to Figs.l9landllOI 



giving rise to a shock surface with an axis ratio of about 3:1. 
At t - 49.46 ms, magnetic (hoop) stresses in the pinched 
toroidal field have accelerated the gas further yielding a so- 
called "nose cone" (well known from simulations of magne- 
tized jets) that is visible in the outermost part of the jet-like 
outflow (Fig. in}. 

The GW ampHtude of model A3B3G5-D3M13 (Fig.[3, 
whose weak-field counterpart A3B3G5-D3M10 emits a GW 
signal of type 111 (Fig.|2j, is modified by the strong initial field. 
Relative to the weak-field model, the pre-bounce maximum is 
enhanced, and the size of the (negative) peak at bounce is re- 
duced by the same amount, whereas the amplitude of the first 
(positive) post-bounce peak remains nearly unchanged. As the 
core loses angular momentum, the size of the hydrodynamic 
and gravitational contributions to the GW amplitude decreases. 
The post-bounce oscillations of the signal, which are mainly 
due to the central core, are superimposed on a nearly con- 
stant positive amplitude of half the size of the bounce ampli- 
tude indicating the presence of the collimated bipolar outflow 
(Fig.ll3>. The evolution of the magnetic contributions of mod- 
els A1B3G3-D3M13 and A3B3G5-D3M13 are similar. In both 
models, a phase of positive amplitude concurrent with the most 
efficient braking of the core's rotation is followed by a decrease 
to large negative values (Fig.^^ lower right panel). In the latter 
phase the magnetic amplitude, which is mostly produced in the 
central core, provides the main contribution to the total signal. 

Increasing the initial magnetic field strength to ~ lO'^G in 
models which exhibit a type-1 GW signal in the non-magnetic 
case, the size of the large (negative) bounce amplitude de- 



creases slightly. This also holds for even stronger magnetic 
fields in case of models that are influenced by rotation to a 
higher degree such as A3B2G4-D3Mm and A3B3G4-D3Mm, 
respectively. However, for rigidly and moderately fast rotat- 
ing models (AlBlG3-D3Mm and AlB3G3-D3Mm) the size 
of the bounce signal grows again strongly when the initial field 
strength reaches lO'^G. 

The ratio of the GW amplitudes of the positive peaks im- 
mediately prior to and immediately after the signal minimum 
at bounce can both increase and decrease by the effects of a 
strong magnetic field. In type-I models with a small influence 
of the rotation on the dynamics, a stronger initial field increases 
the ratio, and for sufficiently fast rotating models (A3B2G4- 
D3Mm) it decreases. Two models with a type-I GW signal 
are close to the transition from a bounce caused by pressure 
or centrifugal forces (A2B4G4-DdMm and A3B3G4-DdMm), 
i.e. their cores are only partially stabilized by the stiffening of 
the EOS at nuclear matter density. For both models the am- 
plitude ratio of the pre- and post-bounce peaks increases with 
increasing magnetic field strength, while it decreases for most 
type-Ill models. 

4.2.2. Models bouncing due to centrifugal forces 

The weak-field models of series A3B3G3-D3Mm belong to 
the transition class between standard-type and centrifugally 
bouncing cores. They bounce mostly due to centrifugal forces 
and exhibit large-scale pulsations after bounce, but their max- 
imum density exceeds nuclear matter density during bounce. 




Fig. 14. The dynamical evolution and the GW signal of model A3B3G5-D3M13. The upper panels display the temporal evolution 
of the rotational energy parameter j0,ot (left), the maximum density p,^^^ (middle), and the magnetic energy parameters (right; 
jSmag (solid line), yS^ (dashed line), and ySmag - yS0 (dotted line)), respectively. The GW signal of the model is shown in the lower 
panels: total amplitude (solid line, left), -Af^.j^^^ (solid line, middle), A^q.^^.^^ (dashed line, middle), and A^q.^^^^ (right). 



Both the period and the damping of the pulsations are signif- 
icantly larger than for purely centrifugally bouncing type-II 
models (Fig.ll5>. Their GW signal consists of a pronounced 
peak of negative amplitude at bounce, and subsequent rela- 
tively long-period oscillations. 

For the strong-field variants of this model series, e.g. model 
A3B3G3-D3M12 (Fig.lB), the main evolutionary effect of 
the magnetic field is the extraction of rotational energy after 
bounce that proceeds on time scales much longer than the dy- 
namic time scale of the inner core. Model A3B3G3-D3M12 
reaches a maximum rotation parameter yS™" = 0.16 at bounce, 
which is nearly the same value as in the non-magnetic case. At 
the time of the second density maximum ySror is only slightly 
smaller than in the weak-field model A3B3G3-D3M10, but at 
the next (3rd) maximum the rotational energy has decreased by 
about 15 %. Unlike in the weak-field model no fourth density 
maximum occurs. Instead, /3roi decreases monotonically. The 
rotation rate of model A3B3G3-D3M12 is sufficiently large 
for a sufliciently long time to allow the core to exhibit sev- 
eral centrifugal pulsations before it settles down in a pressure 
supported final equilibrium state. During these pulsations the 
time-averaged value of p^ax increases, and densities signifi- 
cantly larger than the bounce density ph = 2.41 ■ 10''* cms"' 
are reached. 



The GW signal of the strong-field model A3B3G3-D3M12 
is similar to that of the weak-field model A3B3G3-D3M10 
during the first » 10 ms after bounce when the cores of both 
models still undergo large-scale pulsations. Later when the ro- 
tational energy has decreased sufficiently and the pulsations 
begin to fade away, the signals begin to differ slightly. The 
maximum density now increases monotonically. Rotation is no 
longer important for the core's stabilization, as it is supported 
by pressure forces. Thus, like in the case of a standard-type 
core (e.g. AlB3G3-DdMm), the rapid oscillations of the GW 
amplitude for f ^ 60 ms occur on the core's dynamic time 
scale. Superimposed to the oscillations is a positive mean am- 
plitude ^ 100cm, which results from the enhanced as- 
phericity of the shock wave as compared to that of the weak- 
field model. 

The most strongly magnetized model of series A3B3G3- 
D3Mm, namely model A3B3G3-D3M13, behaves very dif- 
ferently from its weak-field variants as e.g. model A3B3G3- 
D3M10 (Fig.ll5>. Unlike the latter model it suffers a single 
bounce at pi, = 2.69 ■ 10''* cms"' (a value about 10% larger 
than for the weak-field models) with a subsequent contrac- 
tion phase caused by the decrease of the core's rotational en- 
ergy. After bounce the maximum density reaches a value of 
= 3.3 ■ 10"^ cms-'. Unlike model A3B3G3-D3M12, it 
does not exhibit any large-scale pulsations. Instead, it rapidly 
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Fig. 15. The evolution of the rotational energy parameter jSror (left panels), of the maximum density (middle panels) and the GW 
signal (right panels) of models A3B3G3-D3M10 (top panels), and A3B3G3-D3M12 (bottom panels), respectively. 



approaches a high density state without an intermediate phase 
where the density in the entire core is less that nuclear mat- 
ter density, as it is the case for model A3B3G3-D3M12. The 
rotational energy (/Smax = 0.135) is less than for models with 
weaker initial fields (e.g. yS^ax - 0.16 for model A3B3G3- 
D3M10), and decreases by 90 % within 15 ms (Fig.[l6j. 

The collapse and post-bounce dynamics of model 
A3B3G3-D3M13 is similar to that of strong-field single- 
bounce models, as e.g. model A1B3G3-D3M13. Many features 
of the latter model are also present in model A3B3G3-D3M13: 
the bulb-like shock surface, the two highly magnetized post- 
shock regions (one near the equator and the other one along the 
polar axis), and a region of retrograde rotation near the equato- 
rial plane at the edge of the inner core. The ratio of magnetic 
and gas pressure is largest at the axis at r = 175 km well behind 
the shock located at rshock ~ 200 km. 

The GW signal of model A3B3G3-D3M13 is quahtatively 
similar to that of its less strongly magnetized variant A3B3G3- 
D3M12, but the features caused by the magnetic field are 
more pronounced (Fig.ll6>. The peak amplitude at bounce 
A^^ = -1 128 cm is, contrary to model A3B3G3-D3M12, con- 
siderably less negative than in the weak-field model A3B3G3- 
D3M10 (A^^ = -1401 cm). Immediately after bounce the GW 
signal shows certain similarities with that of model A1B3G3- 
D3M13 (Fig.|8}: a large positive peak is followed by a series 
of oscillations occurring on the local dynamic time scale su- 
perimposed on an increase of the mean amplitude to a value of 



~ 600 cm within f « 5 ms, which is the time it takes to 
turn the initially roughly spherical shock into a relatively wide 
bipolar outflow directed along the rotation axis. 

Similar to cores bouncing due to pressure forces, the extrac- 
tion of rotational energy from cores bouncing due to centrifu- 
gal forces proceeds qualitatively differently for models with an 
initial magnetic field strengths of 10'^ G and 10'^ G, respec- 
tively. For the weaker field the extraction process relies on an 
instability-driven angular momentum transport, which is not 
required for the stronger field models. 

Unlike the other models, which collapse to higher den- 
sities when angular momentum transport by the magnetic 
field becomes dynamically important, model A2B4G1-D3M13 
bounces at a lower density (pb = 6.84 ■ 10'^ cms"' compared 
to pb = 1.14 ■ lO'-' cms"' in the weak-field case). At bounce 
its rotational energy (/Jmax = 0.093) is less than that of model 
A2B4G1-D3M11 (J3^^^ = 0.118), which has an initially 100 
times weaker magnetic field. After bounce the rotational en- 
ergy continues to decrease, and the core eventually suff'ers a 
second collapse that is stopped at p^ax - 2.5 ■ 10'^ cms"' 
(Fig.ll7>. During the second collapse /3i-ot rises temporarily 
slightly again, but angular momentum redistribution prevents 
the shock formed by centrifugal forces (at a similar density as 
the one resulting from the first bounce) from stopping the in- 
fall. Thus, a large fraction of the inner core continues to fall 
towards the center, until the collapse is halted by the stiffen- 
ing of the equation of state, i.e. by a bounce due to pressure 




Fig. 16. The dynamical evolution and the GW signal of model A3B3G3-D3M13. The upper panels display the temporal evolution 
of the rotational energy parameter j0,ot (left), the maximum density p,^^^ (middle), and the magnetic energy parameters (right; 
jSmag (solid line), (dashed line), and jSmag - y60 (dotted line)), respectively. The GW signal of the model is shown in the lower 
panels: total amplitude (solid line, left; the dashed line gives the contribution of the layers with r < 55.6 km), -Ajo.fj^j (solid line, 
middle), A^^.^^^^ (dashed line, middle), and Ajo-n^^, (right). 



forces. A shock forms at the edge of the still very extended in- 
ner core (rshock - 70 km compared to rshock ~ 10 km in the 
previously discussed models) at f = 108 ms at the beginning of 
the "plateau" in the evolution of p^ax (see upper middle panel 
of Fig. ll7l l. and it becomes slightly prolate. About 15 ms later 
(f - 123 ms) this shock has reached a radius of rshock ~ 200 km 
(Fig.[T8j. The magnetic field of model A2B4G1-D3M13 is 
less amplified and deformed than in the deeper collapsing 
single-bounce models such as e.g. model A1B3G3-D3M13. 
The largest magnetic to gas pressure ratios are reached well 
behind the shock at about 0.75 rshock ~ 400 km. 

The GW signal does not resemble any of the types used 
to classify non-magnetic models. Hence, we suggest to intro- 
duce a new type-IV GW signal (Fig. ll7l i which it is weaker 
at bounce (A^q = -440 cm) than in the non-magnetic case 
(Ajg = -610 cm), but the amplitude remains large for a longer 
time period due to the longer phase of maximum contrac- 
tion. Until the second collapse the signal varies on time scales 
(Af Si 5 ms) that are short compared to the pulsation periods of 
the non-magnetic case, but comparable to the sound crossing 
time of the relatively extended inner core. After t ^ 140 ms 
the oscillation frequency and the amplitude of the GW signal 
strongly increases as the second collapse begins. During the en- 



tire post-bounce evolution the large negative gravitational and 
hydrodynamic (dominated by the centrifugal amplitude) contri- 
butions to the GW signal, and the large positive magnetic con- 
tribution add up to a total GW amplitude which is very small 
compared to the individual contributions. 

Model A4B5G5-D3M13 bounces at a slighdy higher max- 
imum density (pb = 2.1 ■ 10'"* cm s"') than its weak-field vari- 
ant A4B5G5-D3M10 (pb = 1-97 ■ lO''* cms"'; Fig.O. The 
maximum rotational energy at bounce is barely affected by the 
presence of the strong magnetic field, but it changes the den- 
sity structure of the core completely in the subsequent evolu- 
tion. Models A4B5G5-D3Mm (m < 12) maintain a toroidal 
density structure (rtor ~ 40 km) during their entire evolution 
(Sect.l4.1>. but in model A4B5G5-D3M13 the density struc- 
ture is changed by magnetic forces. Up to a few milliseconds 
after bounce (t $ 32 ms), the core of model A4B5G5-D3M13 
evolves only marginally differently than its weak-field variants. 
Later, however, magnetic stresses transfer angular momentum 
from the still toroidal density distribution into the surrounding 
gas, whereby the torus transforms into a flattened configura- 
tion, which has its density maximum only slightly off-center 
(Fig. l20t . During this phase, the core develops a retrograde ro- 
tating region near the equatorial plane. The core's transforma- 
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Fig. 17. The dynamical evolution and the GW signal of model A2B4G1-D3M13. The upper panels display the temporal evolution 
of the rotational energy parameter j0,ot (left), the maximum density p,^^^ (middle), and the magnetic energy parameters (right; 
jSmag (solid line), yS^ (dashed line), and ySmag - yS0 (dotted line)), respectively. The GW signal of the model is shown in the lower 
panels: total amplitude (solid line, left), -Af^.j^^^ (solid line, middle), A^q.^^.^^ (dashed line, middle), and A^q.^^^^ (right). 



tion process is associated with a large increase of its maximum 
density comparable with that occurring during the first collapse 
(Fig.d. 

During the second collapse, both the magnetic and the rota- 
tional energy of the core increase. The magnetic energy (j8,nag = 
0.11) exceeds that of the primary collapse (ySmaa = 0.078), 
and the rotational energy exceeds slightly the limit for the dy- 
namic instability of MacLaurin ellipsoids (jSrot ~ 0.275), but 
only for about 1 ms. We note that the rotational energy reaches 
its maximum value well before the density does, which is dif- 
ferent from all other models. Later on, the rotation parameter 
rapidly decreases well below the corresponding value of the 
non-magnetic case. 

The shock has a very prolate shape at core bounce, its axis 
ratio being ^ 2 : 1 . At large latitudes, the ratio of magnetic and 
gas pressure is largest well behind the shock (r < 0.75 /"shock ~ 
110 km; Fig.l21t. Due to the rather extreme rotation of the 
model, the high latitudes of the core are strongly rarefied, and 
the shock propagates into a very thin medium near the axis, its 
structure remaining unchanged. 

The GW amplitude of the weakly magnetized non- 
magnetic) models A4B5G5-D3Mm (m < 12) is characterized 
by a very large negative amplitude at bounce (A^q ~ -4100 cm; 
Fig.[l9land Table lD.ll . The signal is positive for several mil- 
liseconds after bounce due to the aspherical shock wave, and 



approaches zero rapidly after the violent re-expansion of the 
core. The bounce amplitude is significantly lowered (A^^ ~ 
-3470 cm) in the strongly magnetized model A4B5G5-D3M13 
(Fig.^]and Table lD.ll . immediately after bounce the GW am- 
plitude is strongly positive and eventually becomes even more 
positive due to the growing asphericity of the shock wave. The 
secondary collapse of this models shows up in the GW ampli- 
tude in the form of a local minimum. Afterwards the amplitude 
rises to values A^^ > 2000 cm and shows superimposed very 
weak oscillations with a period of slightly less than 1 ms. In 
this model, the magnetic and the combined hydrodynamic plus 
gravitational amplitude contributions show a clear phase shift 
resulting from the opposing actions of hydrodynamic (mainly 
centrifugal) and magnetic forces. 



In all our models bouncing due to centrifugal forces the 
peak value of the GW signal at bounce decreases with increas- 
ing magnetic field strength, and for most models the ratio of 
the amplitudes of the post-bounce to the pre-bounce signal 
peaks decreases. Note that the latter statement must be consid- 
ered carefully since - as discussed above - the GW signal may 
be subject to large modifications for models with very strong 
fields. 
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Fig. 19. Evolution of the rotational energy parameter jSrot (left) and the maximum density (middle) of models A4B5G5-D3M10 
(upper row) and A4B5G5-D3M13 (middle row), respectively. The total GW amplitude of model A4B5G5-D3M10 is shown in the 
upper right panel, and the evolution of the magnetic energy parameters, ySmag (solid), yS^ (dashed), and yS^ag -yS^ in the middle right 
panel, respectively. The panels in the bottom row show different contributions to the GW amplitude of model A4B5G5-D3M13: 
total amplitude ( left), -A^l^ ^ (solid line, middle), A^q^ (dashed line, middle), and A^q. (right). 



4.3. Magneto-rotational instability 

The magneto-rotational instability (MRl) is a shear instability 
occurring in differentially rotating magnetized plasma, which 
generates turbulence that am plifies the magnetic fie l d and 
transf ers angular momentum jBalbus & Hawlevlll99l[ 1 19921 
Il998h . If effects due to buoyancy are neglected, a linear sta- 
bility analysis shows that the condition for the instability of a 
mode with wavenumber k is 



dQ2 
dlncT 



+ (cAkr < , 



(20) 



where Q is the angular velocity, m the distance from the rota- 
tion axis, and Ca = b/ yfp the Alfven velocity. When the mag- 
netic field is very small (i.e. the Alfven velocity is very small 
compared to both the local sound speed as well as the local 
rotation velocity) and/or the wavelength is very long, (caA:) is 
negligible, i.e. the MRI occurs simply when the angular veloc- 
ity gradient is negative llBalbus & Hawlevfl99llll99i) : 



dQ2 
din tir 



< 0. 



(21) 
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Fig. 20. Comparison of two snapshots of model A4B5G5-D3M13. The innermost region of the core is displayed at f = 30.98 ms 
(at core bounce; left), and at f = 42.54 ms (right). The logarithm of the density (color-coded), the velocity (arrow), and the 
magnetic field lines are displayed, respectively. 




Fig. 21. Model A4B5G5-D3M13 at f = 30.98 ms (core bounce). The flow field (vectors) and the magnetic field lines are 
displayed together with the gas pressure (color-coded, left) and the magnetic pressure (color-coded, right), respectively. 
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Fig. 12. The hydrodynamic plus gravitational parts of the GW 
amplitude of model A1B3G3-D3M13 at t = 68.5 ms. Positive 
contributions to the integrand in the formula for the hydrody- 
namic plus gravitational amplitude are displayed by colors. The 
bipolar outflow shows up in the large positive parts near the 
axis. 



The time scale of the fastest growing unstable mode is 



■ An 



dO 



dlncT 



(22) 



which depends neither on the strength nor on the geometry of 
the magnetic field. 

General theoretical considerations and non-linear simula- 
tions show that the magnetic energy achievable by the MRI 
field amplification process is of the order of the rotational en- 
ergy, which is comparable to the saturation field expected from 
the process of winding-up field lines by differential rotation 
jAkivama et alJ2003h . However, in axisymmetric (non-linear) 
hydrodynamic simulations of accretion disks the field built up 
by the MRI was found to decay with a rate depending on the 
dissipation properties o f the numerical scheme, and par ticularly 
on the grid resolution jBalbus & Hawlevl ll99l'.'l998'). This is 
due to Cowling's anti-dynamo theorem (Shercliff 1965) ac- 
cording to which an axisymmetric dynamo cannot work in a 
dissipative system, i.e. three-dimensional simulations are re- 
quired. 

In a star a convectively stable stratification will tend to sta- 
bilize the MRI, while convective instability will strengthen the 
MRI. Including the effects of buoyancy a mode with wave vec- 
tor k is unstable with respect to the MRI in a weakly magne- 
tized rotating system, if the instability criterion 
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Fig. 13. Snapshot of model A3B3G5-D3M13 at f = 49.46 ms 
showing the velocities of inflows (blue colors) and outflows 
(red colors). 
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Fig. 18. Model A2B4G1-D3M13 at t = 123.39 ms when the 
expansion of the core is maximum just prior to its second col- 
lapse. Besides the magnetic pressure (color shaded), the flow 
field (vectors), and the magnetic field lines are shown. 



is fulfilled (iBalbus & Hawlevlll99l[rBalbuslll995l) . where N is 
the Brunt- Vaisala or buoyancy frequency. The wave number of 
the fastest growing mode is given by 



(23) (caA:)Lx = 



2\2 



1 - 



16Q4 



(24) 
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where = AD? + dQ^/dln m is the epicycHc frequency. This 
mode grows exponentially with the time scale 



AO. 



in dlntiT 



(25) 



which is a generalization of Ea.E^l 

For decreasing wavenumber k (provided that kc^ < 
(kcA)m-dx) the maximum growth rate {- In/Tmax) decreases ap- 
proaching the value 

<g = -A^'-^' (26) 

in the limit of very long wavelengths (k 0) independently 
of the Alfven velocity. If the MRI occurs in a buoyantly sta- 
ble stratification (where A^^ > 0), long modes will grow very 
slowly compared to the maximum growth rate. In the oppo- 
site case modes with a very small wave number grow very 
fast, and the growth rate of long modes will be of the order of 
the maximum growth rate. In the following discussion we use 
the names "magneto-shear" and "magneto-convective" limit 
to discriminate the two limiting cases without implying any 
deeper physical significance, because it is a matter of ongoing 
debate whether there exists a fundamental physical distinction 
between the two cases in the context of angular momentum 
transport in accretio n disks (see, e.g. (Ba lbus & Hawlevll2002t 
iNaravan et al .P2002')). 

For a rotating magnetized configuration with axial symme- 
try the growth time scale of the maximum growing MRI mode 
in the presence of entropy gradients is given in spherical coor- 
dinates (r, ff) bV iAMy ama et al, (2003) 



2 n 



(jf -2j]+ 1)Q^ + ^— ^ + J] 



16Q2 



da 



2 \^ 



din r 



2 

-1/2 



din r 



where 



Tf = (1 -siniey , 

f = sin^e (1 -sin26') 



(27) 



(28) 
(29) 



The additional subscript "en" refers to the entropy gradient 
included into Eq.|^via the Brunt-Vaisala frequency. In the 
equatorial plane, where the polar angle - n/l, the time scale 
given by Eq.|^is equivalent to that obtained from Eg .1251 

In order to see the MRI in our simulations, we must have 
sufficient spatial resolution to resolve at least the modes with 
the longest wavelengths which grow fastest. However, as these 
modes have the smallest wavenumbers, the Alfven velocity and 
thus the magnetic field (parallel to the direction of k) must 
be sufficiently large for the product Ca ■ k to be in the range 
of maximum growth. This is unproblematic in the magneto- 
convective case due to ffie large growth rate in ffie limit c^k — > 
0. In the magneto-shear limit only modes near (cAk)^^^ grow 
fast, which are, however, difficult to resolve. Further numeri- 
cal complications arise from the fact that we have to identify 
the MRI in a very inhomogeneous and highly dynamical back- 
ground flow. 



Apart from the direct solution of the MHD equations, there 
are alternative ways to investigate the MRI and the turbulence 
driven by the instability, such as the inclusion of well-suited 
models for turbulent transport coefficients into the momentum 
and energy equations. Various closu re models for magnet oro- 
tational turbulence exist. The ones bv lOgilvi? (2001. 2003) and 
IWiUiams (.2005i..2004) which make use of ffie analogy of MHD 
turbulence wiffi viscoelastic flows as observed, e. g. in poly- 
mer fluids in the laboratory, seem particularly interesting. Such 
methods provide a way for including turbulence effects into 
a numerical simulation that cannot treat these effects due to, 
e. g. resolution or symmetry constraints. However, in the simu- 
lations we report here, we did not consider any of these models. 
Instead, we focussed on the possibility of directly simulating 
the development of the instabilities. 

In most of our models we find extended MRI unstable re- 
gions at various epochs. During collapse the cores are unsta- 
ble in the magneto-shear limit. However, the growth times are 
significantly larger than the collapse time scale: even for the 
fastest growing mode they are of the order of a second (in the 
initial models). As the rotational energy and the degree of dif- 
ferential rotation increase during collapse, the growffi times be- 
come smaller, but remains larger than ffie time until bounce, 
i.e. the dynamical background evolves faster than any unstable 
mode. Hence, we do not observe the growth of the MRI even 
for models where the initial magnetic field is sufficiently strong 
for magneto-shear modes to be numerically resolved. 

During and shortly after bounce the cores possess a con- 
vectively stable stratification, i.e. the MRI is of magneto-shear 
character. In model A1B3G3-D3M10 the immediate post- 
shock region is unstable shortly after shock formation due to a 
large negative angular velocity gradient (in ijT-direction). The 
growth times are in the sub-millisecond to millisecond range, 
i.e. comparable to the dynamic time scales. The fastest grow- 
ing modes correspond to spatial scales of less ffian 100 m, 
which are significantly smaller than the (local) grid resolu- 
tion of ~ 300 m. For models wiffi a stronger magnetic field 
(e.g. A1B3G3-D3MI2) ffie spatial scales of ffie fastest growing 
modes are larger due to the larger Alfven velocity, and they can 
be resolved. However, they cannot be discriminated from the 
dynamic background flow at the time of bounce which is dom- 
inated by numerous pressure waves propagating through the 
inner core and launching the shock wave. Moreover, obscur- 
ing the action of the MRI, the strong shear flow suiTounding 
the "surface" of the inner core produces a strong toroidal mag- 
netic field component. The positive entropy gradient behind the 
shock wave does not allow the MRI to grow in most of the 
post-shock region, while large parts of the inner core are MRI 
unstable throughout ffie entire post-bounce evolution due to the 
very flat entropy profile of the core. However, as the instability 
is of magneto-shear nature, one encounters the previously dis- 
cussed spatial resolution problem unless the magnetic field is 
already quite strong. The maximum growth rates are moderate 
(fgi-ow ^ 10 ms). 

Well after bounce when the shock wave begins to weaken 
and thus (locally) negative entropy gradients develop, the 
MRI growth rates become large even for long wavelengths 
(magneto-convective limit). Thus, we observe ffie emergence 
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of the corresponding modes, and a growth of the vorticity oj - 
V X V of the flow. However, this occurs late after bounce, when 
our assumptions concerning the microphysics do no longer 
hold. For model A1B3G3-D3M10, we can study the onset of 
magneto-convective modes. In regions where the growth times 
of the MRI modes (in the limit ^ — > 0) are less than 10 ms 
(Fig. 123, we find a rapid increase of the magnetic energy, the 
magnetic energy density, and of the specific magnetic energy. 
Their temporal variation can be described by exponential laws 
in the interval 57 ms < f < 96 ms. In phases when the field 
is amplified by differential rotation only, the toroidal field in- 
creases but the poloidal field and the corresponding magnetic 
field energy stay approximately constant (Fig. l22t . However, 
when the MRI sets in, the poloidal component grows too due to 
the stretching of the poloidal field lines by meridional motions. 
Because of this behavior the two field amplification mecha- 
nisms can be distinguished. 

For initially very strong magnetic fields (^ 10'^ G), the spa- 
tial scales of interest for the MRI can easily be resolved even 
in the magneto-shear limit, and the MRI criterion is fulfilled 
the growth rates being similar to those of the weak-field case. 
The proto-typical model A1B3G3-D3M12 shows the creation 
of a large amount of vorticity both in the post-shock gas and in 
the MRI unstable regions of the inner core (Fig.Q. Its vortic- 
ity is considerably larger than in the corresponding weak field 
model A1B3G3-D3M10, where we cannot resolve the growth 
of the unstable modes due to their very small spatial scales. 
The flow is organized in sheet-like structures which are typi- 
cally a few kilometers wide in tu-direction and ~ 10 km long 
in z-direction. Another prominent flow structure forms near the 
equatorial plane in the outer layers of the inner core, where a 
violent meridional flow leads to the largest amplification of the 
poloidal field and to angular momentum transport out of the 
inner core. The latter, which implies that the loss of the rota- 
tional support of the inner core is caused by the MRI, holds 
for all models where we can resolve the unstable modes, and 
where the field is not too strong initially. For the extreme mod- 
els AaBbGg-D3M13 the stress tensor is initially already large 
enough for significant angular momentum transport to occur, 
i.e. only a very small amount of field amplification is required. 

Initially weak magnetic fields should lead to a similar dy- 
namical behavior, as the MRI does not depend on the initial 
field strength. The modes should grow on similar time scales, 
but smaller spatial scales. In practice, however, structures re- 
solvable on our grid are much larger than the wavelengths 
of the fastest growing modes, i.e. they have relatively small 
growth rates. Hence, the development of the MRI is artificially 
delayed. In addition, the smaller scales of the more unstable 
modes may affect the global dynamics. Indeed, if we simulate 
the evolution of model A1B3G3-D3M1 1 for a longer time, the 
poloidal field energy starts to rise after f ^ 65 ms. This rise is 
accompanied by the formation of a flow pattern consisting of 
sheets of large vorticity, and the extraction of rotational energy 
from the core, quite similar to our findings for model A1B3G3- 
D3M12, though less pronounced. Thus, the distinction between 
weak- and strong-field models may be irrelevant, because even 
small initial fields give rise to the same kind of instability and 
can become "strong" according to our criterion very rapidly. 



The poloidal field strength oscillates similarly as the den- 
sity or the rotational energy during the post-bounce phase for 
the models bouncing due to centrifugal forces. The field is sub- 
ject to compression and expansion, but apart from the oscilla- 
tions it does not grow unless an MRJ-like mechanism sets in. 

The magneto-convective motions found during the later 
stages of weak field models are much less pronounced in the 
cores of strong field ones. There are two reasons for this behav- 
ior: Firstly, due to the extraction of angular momentum from 
the core one of the driving forces of the instability is reduced, 
and secondly a strong magnetic field is known to suppress con- 
vection. The latter also holds for the MRI which is hampered by 
too strong fields, as the seed field cannot be amplified beyond 
the MRI saturation value. This can be seen from Eq. M3\ : For 
large Alfven velocities only modes with very long wavelengths 
(or small wavenumbers) are unstable. If the critical wavelength 
exceeds the intrinsic scales of the problem, e.g. the size of the 
inner core or the size of potentially unstable regions, the un- 
stable longer modes cannot develop, and the modes of smaller 
wavelengths are stable. 

4.4. Influence of the magnetic-field configuration 

The dynamics and gravitational wave signature of the mod- 
els A3B3G5-DdM12 differing only in the location of the 
field-generating current loop (i.e. in their Dd parameter; see 
Sect. 12. 3.^ are qualitatively quite similar, except for model 
A3B3G5-D0M12 whose current loop is located at infinity (see 
the discussion below). The time scale for the slow-down of 
the core exhibits some variation, but otherwise the evolution 
and the GW signal agree qualitatively quite well with those 
of the corresponding reference models A3B3G5-D3Mm (see 
above). This may be surprising because the initial total mag- 
netic energies of models having uniform fields may be more 
than 100 times greater than that of current-loop models (table 
13 • However, the outer parts of the very extended field struc- 
ture of the more strongly magnetized models, where most of 
the additional magnetic energy is stored, are dynamically less 
important than the central regions, and in particular the region 
around the field-generating current loop where the magnetic 
field is strongest. 

For the initially homogeneously magnetized model 
A3B3G5-D0M12 a less efficient braking of the core's rotation 
is observed than for the corresponding current-loop models. 
Instead of a contracting secularly, the core expands to sub- 
nuclear densities on a comparable time scale. The rotation en- 
ergy of the immediate post-bounce core of model A3B3G5- 
D0M12 is sUghtly smafler (J3,ot ~ 0.076 at f » 32 ms) than 
for the corresponding non-magnetic core (fii-ot ~ 0.084 at 
f s; 32 ms), and it decreases sightly during the subsequent 
evolution. Although the magnetic field gets amplified a little 
during collapse, ySmag decreases until bounce because most effi- 
cient amplification takes place in the central regions of the core 
whereas the total magnetic energy is dominated by the outer 
regions where the field stays roughly constant during collapse. 
Afterwards, both E^ng and /Smug show a comparable growth to 
the current-loop models, and the ratio of toroidal field energy 



M. Obergaulinger et al.: Magneto-rotational core collapse 



25 



log p [cgsj 

il.iL 



i:.z3 



13.34 



14.^ 






Fig. 23. The post-shock region of models A3B3G5-D0M12 (left quadrants) and A3B3G5-D2M12 (right quadrants) at t ^ 
33.6 ms about 3 ms after core bounce. Density and magnetic pressure (color-coded) are displayed in the top and bottom quadrants, 
respectively. Additionally, the flow field (vectors) and the poloidal field lines of the magnetic field are displayed. 



and total field energy grows steadily, but slowly. The magnetic 
field of the uniformly magnetized core has a simpler struc- 
ture at late times, while the current-loop fields become more 
and more twisted. The gravitational wave amplitude of model 
A3B3G5-D0M12 is strongly dominated by the fluid contribu- 
tion in spite of its large initial magnetic energy, and resem- 
bles that of the corresponding current-loop cores in many re- 
spects. The small magnetic contribution to the signal never ex- 
ceeds that of model A3B3G5-D3M12, although this model has 
a 15 times smaller magnetic energy (but the same central field 
strength). 

Model A3B3G5-D2M12 has an initial magnetic field en- 
ergy which is a factor of 150 greater than that of model 
A3B3G5-D0M12 (table|3. Nevertheless, at f a; 39 ms, its mag- 



netic field is almost three times as energetic as that of model 
A3B3G5-D0M12 (/3mag = 0.0075 compared to y3n,ag = 0.0028) 
due to the different amount of field amplification occurring in 
both cases. At bounce both models have developed a region 
of high magnetization near the edge of the inner core close to 
the rotation axis. In model A3B3G5-D2M12 this region also 
extends towards larger radii along the equatorial plane form- 
ing a highly magnetized "sheet" which is a relic of the initially 
highly magnetized region around the field generating current- 
loop. The shape of the shock waves is similar in both models, 
but the regions of large magnetic pressure are more extended 
in the current-loop model filling the whole post-shock region, 
whereas they are more concentrated towards the very center of 
the core for the uniform field model (Fig.|23}- 
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Fig. 24. T[ie evolution of the maximum density (top panels) and the GW amplitude (bottom panels) of models A3B3G3-D0M13 
(left panels), A3B3G3-D3M13 (middle panels), and A3B3G3-D1M13 (right panels), respectively. The GW signal of model 
A3B3G3-D0M13 shows a constant offset from zero, which is subtracted here. 



Similar trends as for model series A3B3G5-DdM12 can 
be observed in the case of model series A3B3G3-DdM13 
(Fig.l24>. The dynamics of the current-loop model A3B3G3- 
D1M13 is qualitatively comparable to that of a A3B3G3- 
D3Mm model with a much weaker magnetic field. Hence, its 
GW signal does not differ much from the non-magnetic case 
in spite of its strong initial magnetic field. However, we find 
large differences between the uniform field model A3B3G3- 
D0M13 (Fig.|23 left panels) and the current-loop models, as 
e.g. A3B3G3-D3M13 (Fig.|^ middle panels). The maximum 
density of the former model exhibits large-scale pulsations dur- 
ing which the bounce density is exceeded by 9 %, and p^ax is 
smaller than pnuc for most of the time. The core has a very high 
magnetic energy initially, ySjnag ~ 20 %, which decreases during 
the collapse until it levels off at a; 7 %. This value is still much 
larger than the largest one reached by the current-loop mod- 
els during their whole evolution. The GW signature of model 
A3B3G3-D0M13 reflects the density oscillations of the core, 
and like in the case of model A3B3G3-D3M13, it is consid- 
erably shifted towards positive amplitudes (Fig.l24> due to the 
very aspherical shape of the shock wave. 



5. Summary and conclusions 

Past and present realistic studies of core collapse supernovae 
neglect magnetic fields. Instead, they focus on an elaborate and 



accurate description of the microphysics and (neutrino) trans- 
port physics. Since this is a formidable task requiring extensive 
computational resources, these studies are necessarily limited 
to a very small number of different initial models. Bearing in 
mind our lack of knowledge of the exact conditions in the pre- 
collapse star, we pursued a complimentary approach. Instead of 
simulating the evolution of a few stars in great detail, we per- 
formed a comprehensive parameter study covering the space of 
possible progenitors with a large number of initial models and 
exploring a relatively unknown territory, namely the influence 
of magnetic fields. Consequently, we are restricted to some 
simplified and approximate treatment of the core's physics, and 
employ only simplified progenitor models (i.e. polytropes). 

Our parameter study focuses on the investigation of two 
main questions: 

- How do magnetic fields affect the dynamics of the collapse 
and bounce, and how is their presence reflected in the grav- 
itational wave signal emitted by the core? 

- Which pre-collapse magnetic fields (strength and geome- 
try) and angular momentum (amount and distribution) will 
cause important dynamic effects for core collapse? 

To address these questions we have performed the most com- 
prehensive parameter study of magneto-rotational core col- 
lapse up to now. We simulated the gravitational collapse and 
the immediate post-bounce evolution of a series of differen- 
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Fig. 22. Comparison of the evolution of models A1B3G3- 
D3M12 (upper panels) and A1B3G3-D3M10 (lower panels). 
The left panels display the total magnetic energy (solid lines) 
of the core, and the total (dashed lines) and poloidal field en- 
ergy (dash-dotted lines) of the MRI-unstable regions, which 
are defined according to two different criteria: For the strong- 
field model A1B3G3-D3M12 all zones are considered where 
the growth time of the fastest growing mode is less than 10 ms, 
while for the weak-field model A1B3G3-D3M10 the corre- 
sponding magneto-convective limit is used. The energy density 
Cmag - £mag/V of the poloidal field in the MRI-unstable re- 
gions is shown in the right panels. The poloidal field increases 
rapidly when the MRI grows, which happens in the strong-field 
model (top right) shortly after bounce, and in the weak-field 
model (bottom right) by magneto-convective modes at later 
epochs. 

tially rotating, magnetized, axisymmetric stellar core models 
using a newly de veloped 2D New tonian MHD code based on 
the algorithm of Pen et al. ('2003*) that employs the relaxing 
TVD sch eme of Jin & Xin (1995) and the constraint-transport 
method jEvans &Hawlevl[l98i to keep the magnetic field 
divergence-free. Our code, which has been comprehensively 
tested, incorporates an analytic equation of state describing 
in an a pproximate way th e thermodynamic properties of core 
matter (Janka et al. 1993). We also computed the gravitational 
wave signal produced by the core with the Einstein quadrupole 
formula taking into account both the contributions of the hydro - 
dynamic and gravitational forces ( Monchmever et al. '1991"), 
and of the Lorentz forces (Kotake et al. 2004b) . The initial 
models used in our stu dy ar e identical to thos e used b y 
IZwerger & Mulleii ( 1 19971) and 'Pimmelm eier et all ( l2002albl) . 
The initial magnetic field is purely poloidal. 

Our simulations show significant modifications of the dy- 
namics and the GW signal of strongly magnetized cores com- 
pared to non-magnetic ones. Magnetic fields are amphfied in 



some of our models to very large field strengths which are in the 
magnetar range (|B| ~ 10'^ G). Such strong magnetic fields ef- 
ficiently extract rotational energy from the central core, thereby 
triggering a secular post-bounce contraction phase. They also 
cause the formation of collimated outflows, which give rise to 
distinctive features in the GW signal. Hence, it is in principle 
possible to extract information about the degree to which an 
explosion has a jet-like character from the GW signal. Further 
specific results are: 

1. From the initially purely poloidal magnetic field, large 
toroidal field components are created by the action of the 
differentially rotating core in a so-called Q dynamo, the 
energy of the saturation field being of the order of the rota- 
tional energy of the core. The efficiency of this field wind- 
ing process depends on the angular momentum distribution 
of the progenitor and the amount of differential rotation 
arising during collapse. The latter is determined by the de- 
gree of non-homologous collapse which is sensitively in- 
fluenced by the equation of state and the rotation rate. The 
saturation fields estimated from the rotational energy of our 
collapsed cores are of the order of 10'^ G, which are field 
strengths expected for magnetars, but not for typical neu- 
tron stars. The field amplification process is found to be 
less efficient if some magnetic energy is consumed for the 
acceleration of matter by dynamically important fields. 

2. As MHD instabilities acting on the toroidal field compo- 
nent are suppressed due to the assumed axisymmetry, we 
are unable to simulate the corresponding field amplifica- 
tion processes, and thus may underestimate the saturation 
value and the growth rate of the magnetic field. However, 
if the entire dynamo process of transforming poloidal and 
toroidal fields into one another and their amplification is 
mainly driven by differential rotation, the saturation fields 
and the growth rates are expected not to exceed those found 
in our simulations. We point out that (i) without any sym- 
metry restriction fields of similar strengths may arise when 
one evolves initially less magnetized cores, and (ii) that the 
strong initial fields of our progenitors (which are stronger 
than those expected in actual cores) may compensate for 
a possible additional dynamo action unconsidered in our 
models. The field amplification is fastest for rapidly and 
very diff'erentially rotating models. In models producing a 
type-1 GW signal the magnetic field is amplified more effi- 
ciently than both in the faster collapsing models of type-Ill 
and the centrifugally bouncing type-II models. 

3. The influence of a magnetic field on the dynamics and the 
GW signature of the core depends mainly on the initial 
magnetic field strength, and to a lesser extent on the field 
topology. For the weakest of our initial fields (B ~ 10'" G) 
we do not find any significant impact either on the dynam- 
ics or on the GW signal on the time scales covered by 
our simulations 100 ms after core bounce). However, 
sufficiently strong fields have a large impact on the dy- 
namics of the core and on the gravitational wave signal. 
Compared to the coiTesponding non-magnetized models, 
the core bounce is somewhat delayed for the majority of 
our very strongly magnetized models {B S 10'^ G), and 
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the bounce density is slightly enhanced. Extremely strong 
initial magnetic fields (B <: 10'^ G) efficiently slow down 
the rotation rate of the central core already during collapse, 
and thus reduce its. centrifugal support by a considerable 
amount. Due to the loss of angular momentum from its 
central regions the core eventually begins to contract or 
even enters a second collapse phase after bounce. Hence, 
an initially extremely magnetized and rapidly rotating core, 
which bounces due to centrifugal forces at low densities, 
evolves into a very compact configuration that is almost en- 
tirely pressure supported, its initially toroidal density strat- 
ification changing into a flattened centrally concentrated 
one. 

4. The magnetic stresses of the toroidal field component (cre- 
ated by the differential rotation) can cause a strong pinch 
effect, which gives rise to a mildly relativistic, jet-like out- 
flow along the rotation axis, and consequently to a very pro- 
late shock wave. The aspherical outflow is mainly driven by 
the very high magnetic pressure i n the post-shoc k plasma, 
as first pointed out by Yamada & Sawail ll2004l) . and it is 
observed both in cores that bounce due to the stiffening of 
the EOS and (mainly) due to centrifugal forces. In purely 
hydrodynamic calculations highly aspherical outflows are 
only encountered for very extreme rotators, whereas they 
are quite common for models with an initial magnetic field 
strength of Z? ^ lO'^G. On somewhat longer time scales, 
similar but weaker outflows are observed for models having 
only a tenth of this field strength. This finding might be of 
relevance for the study of Gamma-Ray Bursts. According 
to our results extreme initial conditions are required for 
the formation of jet-like outflows on the time scale of a 
few ten milliseconds. However, our results may also ap- 
ply on longer time scales to the smaller rotation rates and 
the weaker magnetic field strengths of more realistic stellar 
progenitors. 

5. All simulated cores develop MRI unstable regions. 
However, in most cases, a fast growth of the MRI is lim- 
ited to short-wavelength modes which we cannot resolve 
due to insufficient numerical resolution. For initial fields 
^ 10'2 G, the fastest growing modes have sufficiently 
long wavelengths to be resolved numericaUy. A meridional 
flow pattern characterized by large vorticity arises the cir- 
culation pattern being organized in sheet-like structures 
which is also the topology of the poloidal field. In con- 
trast to the weak-field case, where we cannot resolve pos- 
sibly unstable modes and where the poloidal field stays 
roughly constant after bounce, the poloidal field energy of 
strongly magnetized cores grows approximately exponen- 
tially. The field strengths found in our simulations are of 
the order of the saturation field estimates derived for the 
MRI, i.e. ~ lO'^ G at the surface of the inner core. The 
growth of the field comes along with the transport of an- 
gular momentum out of the inner core, and thus a loss of 
rotational support. Consequently, the maximum density of 
the core rises steadily during the post-bounce evolution. 

6. The gravitational wave signature of the core is altered 
significantly in the presence of initially strong {B <; 
10'^ G) magnetic fields compared to the coiTesponding 



non-magnetic case. Since pressure-bouncing type-I and 
type-Ill models are affected only quantitatively by strong 
magnetic fields (in contrast to the centrifugally bouncing 
models where magnetic fields change the collapse dynam- 
ics considerably), the GW amplitude of a strongly mag- 
netized type-I or type-Ill model does not deviate strongly 
from the purely hydrodynamic case. Depending on the ac- 
tual rotational profile and the EOS, the GW bounce ampli- 
tude may be enhanced or weakened by a few 1 %. This 
confirms the findings of lYamada & Sawail (l2Q04l) obtained 
for a much smaller set of models. 

7. For most strongly magnetized models we observe a shift of 
the bounce signal towards more positive amplitudes. A ma- 
jor part of this shift is due to magnetic forces: around the 
time of core bounce the magnetic part of the GW amplitude 
is positive and large indicating that the magnetic stresses 
act to diminish the core's oblateness. This reduces the ab- 
solute value of the pronounced (negative) bounce ampli- 
tude for most type-I and type-II models, and enhances the 
(positive) maximum GW amplitude of the type-Ill models 
reached shortly before or after bounce. In the latter models 
strong magnetic fields also decrease the size of the much 
less pronounced (negative) bounce minimum of the GW 
amplitude. Rigidly and modestly fast rotating type-I models 
behave differently, as for these models the size of the (nega- 
tive) bounce peak is enhanced significantly for initially ex- 
tremely strong magnetic fields (~ lO'^G). After bounce, 
both the gravitational and the hydrodynamic parts of the 
GW amplitude decrease as the rotational flattening of the 
central core diminishes due to the extraction of angular mo- 
mentum. Later in the evolution the magnetic part of the GW 
amplitude, now becoming increasingly negative, dominates 
the hydrodynamic one for cores which are slowed down 
very strongly. As in the case of purely hydrodynamic mod- 
els, the GW signal of a strongly magnetized core bouncing 
due to pressure forces exhibits oscillations with periods that 
are comparable to the local dynamic time scale. 

8. The dynamics of type-II models bouncing due to centrifu- 
gal forces is changed completely by the presence of strong 
magnetic fields. This is reflected in a rather drastic change 
of the GW signal emitted by these cores. As the large-scale 
pulsations of the core fade away during the secular con- 
traction phase, the long-period modes of the GW signal 
become less important. Instead high frequency oscillations, 
similar to the ones emitted by type-I models, dominate the 
GW signal. The coiTesponding periods are given by the lo- 
cal dynamic time scales of the core which is no longer sup- 
ported by rotation, but - at much higher densities - by pres- 
sure forces. This shift in the frequency of the GW ampli- 
tude provides observational evidence of the change of the 
core's dynamics and structure by a strong magnetic field. 
We introduce the new GW signal type IV to refer to this 
phenomenon. 

9. Aspherical outflows and shock waves get imprinted on the 
GW signal the signature being a large, slowly varying posi- 
tive GW amplitude. A few (~ 5) milliseconds after bounce, 
we find amplitudes of the order of ~ 10^ cm which re- 
main at roughly this level until the end of our calculations. 
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As bipolar outflows occur preferentially in strongly mag- 
netized cores, the GW signal may be used to distinguish 
strongly magnetized cores from non-magnetic or weakly 
magnetized ones. 

The applicability of our findings to the supernova prob- 
lem is limited by our simplified treatment of the microphysics 
(approximate equation of state), by the neglect of neutrino 
transport, and by our simplified and probably too extremely 
rotating and magnetized initial models. Nevertheless, after 
having shown that our code is suited for multi-dimensional 
MHD simulations of core collapse, and after having high- 
lighted some general trends of magneto-rotational core col- 
lapse, we are ready for a more sophisticated investigation of 
the phenomenon. In particular, the evolution of our (more re- 
alistic) weak-field models on longer time scales will be the 
topic of further studies, as on longer time scales even initially 
weak magnetic fields can be amplified, e.g. by means of the 
magneto-rotational instability, to a strength relevant both for 
the dynamics and the GW signature of the core. Furthermore, 
the interplay of magnetic fields and various hydrodynamic in- 
stabilities occurring in a supernova explosion is a challeng- 
ing topic to be addressed. In particular, the convective region 
between the proto-neutron star and the stalled hydrodynamic 
shock wave might be an arena where generic MHD eff'ects such 
as the magneto-rotational instability do operate. 

Imposing axial (and equatorial) symmetry is a serious re- 
striction when simulating MHD flows. It prevents us from in- 
vestigating the fate of cores whose magnetic field and rotational 
axis are misaligned as in pulsars, and from simulating the in- 
trinsically 3D magneto-rotational instability to the full extent. 
We think that the first shortcoming will influence our results 
only slightly, as the initial magnetic field is soon completely 
dominated in our models by the field generated by differential 
rotation. Thus, provided there is a poloidal seed field of similar 
strength no significant changes are to be expected. The second 
restriction is probably uncritical for the relatively short time 
scales covered by our simulations, as the dominant amplifica- 
tion process for the magnetic field during these early stages 
is the rapid winding-up of the initial poloidal field by means 
of dififerential rotation. Additional field growth by the MRI or 
other MHD instabilities is probably less important. In any case, 
limited by the present computational resources, a 3D MHD pa- 
rameter study of core collapse is still some years ahead. 

A further limitation arises from our Newtonian approach, 
as we observe outflow velocities in some of our models that 
exceed ~ 30 % of the speed of light and that are still in- 
creasing when we had to stop our simulations. Although be- 
ing only mildly relativistic in our models, a study of the further 
evolution of these outflows may require a special relativistic 
MHD code, if the jets continue to accelerate. The eff'ects of 
general relativistic corrections of the gravitational potential on 
magneto-rotational core collapse arising due to the high com- 
pactness of the collapsed cores have been investigated by us in 
some detail. The results of this complementary study will be 
presented in a separate pubUcation. 
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Appendix A: Relaxing TVD 

In this section we summarize the basics of the relaxing TVD 
method. For fu r ther inform atiori the reader is refe iTed to 
IPen et alJll2nn3h : I.Tin & X in ( 1995\: ItY^ & PenI (120031) . 

A major drawback of MUSCL-type and Riemann solver 
schemes for systems of non-linear conservation laws is the 
need for calculating the eigenvalues a,„, and the eigenvectors of 
the Jacobian of the flux vector (Hebrew indices (N = 1 , . . . , m) 
enumerate the components of the system of m equations) 



du- 



(A.l) 



which usually is a computationally quite demanding step. This 
difficult y can b e overcome by employing a relaxing TVD 
scheme il.Tin & Xi n 1995; Trac & Pen 2003 ) that does not re- 
quire the explicit calculation of the eigenvectors of the system, 
and in some cases of most of the eigenvalues. The numeri- 
cal method we have applied turns out to be both robust and 
accurate, and is well suited for simulations of core collapse. 
This is good news, as the method involves several considerable 
simplifications compared to more elaborate methods based on 
Riemann solvers, which make it computationally very cheap, 
and hence attractive. One replaces the one-dimensional system 
of conservation laws 



d,U^ + VF^ = 



(A.2) 



with the initial condition U^{^, 0) = t/y (^) by the relaxation 
system 



d,V^ + A^VU^ = — (y^ - F^) , 

T 



(A.3) 
(A.4) 



with the initial conditions U^i^, 0) = U^i^) and V^i^, 0) = 
F^(Uq(^)). The constant t is called the relaxation rate. 
The constant matrix A = diag(fli, «,„) satisfies the sub- 
characteristic condition 



- /2 > . 



(A.5) 



Thereby, the non-linear system (IA.2> gets replaced by a sys- 
tem of linear equations that can be solved without the need of 
Riemann solvers. In the zero relaxation limit t ^ oo, one ar- 
rives at the relaxed system 



dtV^ +A^VU^ = 0, 



(A.6) 
(A.7) 
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F^{U^). Introducing right (R) and left (L) moving Appendix C: Gravltational-wave extraction 



2 



the system decouples yielding 



(A.8) 



(A.9) 



These two systems (R, L) of equations have the structure of lin- 
ear advection equations for rightward and leftward directed ad- 
vection, respectively. The constant advection velocities in both 
systems of equations are given by 



(A. 10) 



Being of a quite simple structure and having quite simple char- 
acteristics with constant speeds c*'^, the system (IA.9> can eas- 
ily be solved using a first- or second-order accurate upwind 
scheme. The upwind fluxes at the zone interface >S,+i/2 will be 
computed from variables of zone 2/ for the right-moving vari- 
able f/*^'*, and from zone Xi+i for U^'^, respectively. The value 
of can then be recovered from U^''^'^ by 



Since we do not evolve the consistent set of general-relativistic 
MHD equations, we are not able to extract the GW signal from 
the metric. But, in the post-Newtonian limit, it is still possible 
to extract the lowest-order terms of the GW signal from the 
hydrodynamic variables. For this purpose we have used the re- 
formulation of the quadrupole formula by iMonchmever et all 

dim), 

where temporal derivatives are replaced by spatial ones 
using the continuity and the Euler equations. Thus, the hy- 
drodynamic fluxes and source terms are introduced into the 
quadrupole formula. 

As we use spherical coordinates (r, 6, <p) in our simulations, 
the metric perturbation associated with a gravitational wave is 
most conveniently expanded in terms of pure spin tensor har- 
monics r^,^''™ and T^u'''^ with the coiTesponding "electric" and 



ab 



' ah 



"magnetic" amplitudes A^q and Ajq , respectively: 



1=2 m=-l 



+A 



('-!) 

'('-') 



(C.l) 



Under the assumption of axisymmet ry, the only non-vanishin g 
amplitude is A??, which is given by jMonchmever et alJl99lh : 



(A.ll) aEI 



We compare both local properties of the models such as the 
maximum density of the core, and global ones such as the GW 
signal and the evolution of ySrot- The results of these simulations 
indicate that our standard resolution is adequate to follow the 
evolution of our models. 
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As I.Tin&5cinlll 1995 1) showed, the resulting scheme will be 
TVD, if the sub-characteristic condition ( IA.5> . and the CFL 
condition arising form the advection terms are fulfilled. 



Appendix B: Convergence tests 

To ensure the numerical convergence of our results, we per- 
formed a number of simulations of the same model using dif- 
ferent grid spacings. Based on these simulations, we selected 
a grid of 380 zones in radial and 60 zones in angular direction 
as our standard grid. The grid is logarithmically spaced in ra- 
dius, with a central resolution of (Ar)c ~ 300 m and a relative 
increase of the cell size from zone to the zone of » 1.1 %. The 
angular grid is uniform. 

Results from a sample of the convergence runs are dis- 
played in Fig. lB.ll We show the temporal evolution of the max- 
imum density, the GW amplitude, and p^ot for model A1B3G3- 
D3M13, computed with grids of 



- 280 zones in radius and (A,)c ~ 600 m, and 30 zones in 
angle (dotted lines), 

- standard resolution as described above (solid lines), and 

- 580 zones in radius and (A,)c ~ 150 m, and 90 zones in 
angle (dashed lines). 



where the radiative quadrupole moment M|q is (z - cos ff) de- 



fined as 



M^^ - 



G 32n^l^ 



Jo Jo 



3 \2 2 



(C.3) 



Using the continuity equation, the first time derivative of M|q 
is given by 

G 32n^'^ 



E2 
20 



dz d—pr 

Jo Jo J 



Vr\-Z 



l\-3vgz^|l-z^ 



' Vl5 

The time derivative of Nf^ is computed using the MHD mo- 
mentum equation which allows ones to eliminate ^^(pvr,pvg) 
with the help of the hydro-magnetic and gravitational force 
terms given by 

dr^Frr 1 d sin OF re 



F' 



rdi-cosO) 
+ k2P,+pvl-bl+pvl-bl), 



cote 

— (n- 



1 dsm6Fgg 
r d{- cos 0) 



■bl)- 



(pV,.Vg - brbg) 



--pdg^. 
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(C.5) 



(C.6) 
(C.7) 

(C.8) 



The components of the flux tensor in the above expressions are 
given by 



Fab = pVaVb + P*5ab ' babb- 



(C.9) 
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Fig.B.l. The temporal evolution of the maximum density, the GW signal and the rotational energy parameter ySiot of model 
A1B3G3-D3M13 comparing simulations with grids of different resolution. Results from a run with the standard grid of 380 
logarithmically spaced zones in radius (central resolution of 300 m) and 60 zones uniformly distributed over < < n/l are 
shown by the solid lines, and runs using a coarser grid of 280 zones in radius (central resolution of 600 m) and 30 zones in 
angle and a finer grid of 580 zones (central resolution of 150 m) and 90 zones in angle are plotted with dotted and dashed lines, 
respectively. 



Followin g the dei"ivatiq n of the quadrupole formula by 
( Monchmever et al J 1 99 ll) . one then finds 



G 32n-2 



Jo Jo 



CO 



c^'VB Jo '~Jo '3 
fJ3z^-l) + M2-3z^)- 



(>frez Vr 



-rdM3z^ - 1) + 3de<iz Vl - 
where the components of fij are given by 

fij = pViVj - bibj . 



(C.IO) 



(C.ll) 



The total isotropic pressure f * (sum of the gas pressure and 
the magnetic pressure) cancels out, i.e. only the velocities and 
the magnetic stresses appear in the hydro-magnetic part of the 
amplitude. 

Eg. lC. lOl corresponds to the hydrodynamic an d the Lorentz 
force p arts of the quadrupole formula given bv iKotake et alJ 
ll2004bV who derived the following expression for the GW am- 



plitude Af2 = A^l^^^ + AE2^^_.^,^^^ + A^l^^^: 
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dr[Bgr\3z^-\))rE^+ 
+dr{B^r\3z^-l))rEe- 

—^dg(B^sm6r\3z^ - 1))e, 



(C.14) 



To evaluate these expressions, the current density j and the 
electric field E have to be calculated from the magnetic field 
and the velocity. The last term (Eq. lC.14T i describes the con- 
tribut ion resulting from th e energy density of the magnetic 
field. iKotake et alJ ll2004bl) found that the hydrodynamic and 
the Lorentz force contributions are at least two orders of magni- 
tude larger than the magnetic energy one. Therefore, and due to 
the time derivatives still involved in its calculation, we neglect 
this contribution and consider only the GW amplitude resulting 
from the quadrupole moment of the matter. 

Appendix D: Synopsis of our results 



Tables ID. II through ID. 41 provide an overview of the dynamic 
evolution of the flow and the magnetic field, and about the re- 
sulting gravitational wave signal of all our models. 
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Table D.l. Some characteristic model quantities: the first two columns give the model name and the classification of the GW 
signal (for the corresponding non-magnetized model). Columns 3 and 4 give the time of bounce tb (in milliseconds) and the 
maximum density at bounce pb (in units of lO'"^ cms"'). An exclamation mark behind the density value signifies that the maxi- 
mum density of the model exceeds the bounce density during the later evolution. A^q (column 5) and AjQ.^^g (column 6) are the 
maximum GW amplitude (in cm) and the corresponding magnetic contribution. A^q.^ (Column 7) is a rough mean value of the 
wave amplitude (in cm) at some late epoch; no value is provided when the GW amplitude does not approach a quasi-constant 
asymptotic value. If the absolute value of this amphtude is large, the presence of an aspheric outflow at late epochs can be in- 
ferred. The following columns give the maximum value of the rotational (column 8) and the magnetic beta parameter (column 
9), the time when yS^ag reaches its maximum (column 10), and the corresponding beta of the toroidal field (column 11). If the 
magnetic field is still amplifying at the end of the simulation, an exclamation mark is added behind the table entry, and if the 
magnetic field is decreasing at this time, we give its final value p^^^ in parentheses. 
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A3B3G3-D3M10 


n/T 


49.70 


2.41 


-1400 


0.017 


30 


16 


1.3 ■ 10"'! 


75.0 


1.2 ■ 10" 


-3 


A3B3G3-D3M11 


m 


49.70 


2.41 


-1401 


1.5 


30 


16 


0.18! 


73.4 


0.18 




A3B3G3-D3M12 


ii/i 


49.71 


2.42! 


-1379 


56 


40 


16 


14! 


71.1 


11 




A3B3G3-D1M13 


ii/i 


49.79 


2.39 


-1400 


2.2 


40 


16 


0.49! 


77.1 


0.46 




A3B3G3-D3M13 


ii/i 


51.00 


2.68! 


-1128 


705 
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13.5 


5.3 (2.6) 


52.9 


3.1 
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50.41 
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-1339 
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15 


21 (7.4) 
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-4 


A3B3G4-D3M10 
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0.0021 
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A3B3G4-D3M11 


I 


39.64 
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-894.9 


0.21 


15 


14.6 
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69.7 
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I 


39.65 
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65 
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0.034 
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9.5 
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56.3 
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30.35 
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261.6 


0.47 
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1.0 


45.01 
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1.6! 
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-3 
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5.1 (3.1) 


37.4 


2.5 
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h 


Pb 


.a 
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^20;oo 


omax 
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omax (ohn 
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[ms] 


[10'* -^] 


[cm] 


[cm] 


[cm] 


% 


"% 


[ms] 
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A4B5G5-D3M10 


I/II 
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64.0 
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0.69 
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34.6 
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57.6 
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34.5 


0.60! 


32.1 


0.52 
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2.10! 


-3473 


1225 
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Table D.3. Some characteristic model quantities (name of model given in column 1) at time t (in msec; column 2) when the core 
has reached a quasi-equilibrium state. For models which do not reach a quasi-equilibrium state until the end of the simulation 
(e.g. type-II models with large scale core pulsations) we provide upper (top value) and lower (bottom value) bounds estimated 
from the values at maximum and minimum contraction. Columns 3 and 4 give the surface radius (in km) and the mass 
(in solar masses) of the quasi-equilibrium configuration, respectively. Since it is still surrounded by an (expanding) envelope of 
high density matter, the definition of its surface radius Tc is somewhat uncertain. As the rotation rate 2n/Q. (in msec), where Q 
is the angular velocity averaged over the angle 6, as well as the total magnetic field \b\ and (the absolute value of) its toroidal 
component (both in Gauss) vary strongly near the surface and on short time scales, the corresponding values in columns 5, 6 
and 7 should be used with care. Negative values of the rotation rate signify counter-rotating cores. Finally, in columns 8 and 9 
we give the radii of the shock at the polar axis, r^^, and at the equator, r^j^ (both in cm), respectively. No entry in these columns 
impUes that the shock has already left the computational grid. 



Model 


t 






2n/Q. 


\b\ 


1^*1 








[ms] 


[km] 


[Mo] 


[ms] 


[G] 


[G] 


[km] 


[km] 


A1B1G3-D3M10 


75 


22.5 


0.59 


9.6 


4.9 


10'^ 


4.7 


10'^ 


725 


692 


A1B1G3-D3M11 


75 


22.3 


0.59 


8.2 


5.7 


10'* 


5.6 


10" 


725 


692 


A1B1G3-D3M12 


75 


21.5 


0.59 


11.9 


8.7 


10'* 


4.8 


10'* 


725 


692 


A1B1G3-D3M13 


75 


23.4 


0.61 


-104.9 


8.5 


10'* 


4.6 


10'3 


768 


669 


A1B3G1-D3M10 


116 


49.8 


1.2 


8.0 


1.9 


• 10'2 


1.9 


10'- 


920 


745 




126 


129 


1.2 


52.4 


3.6 


• 10" 


3.6 


10" 






A1B3G1-D3M11 


116 


49.7 


1.2 


8.0 


1.9 


10" 


1.9 


10" 


920 


745 




126 


129 


1.2 


46.7 


3.2 


10'2 


3.2 


10'2 






A1B3G1-D3M12 


116 


49.0 


1.2 


7.9 


2.0 


10'* 


1.9 


10'* 


920 


745 




125 


129 


1.2 


45.5 


3.2 


10" 


3.2 


10'3 






A1B3G1-D3M13 


117 


35.3 


1.1 


6.3 


1.4 


10" 


4.7 


10'* 


813 


568 


A1B3G3-D3M10 


70 


25.5 


0.68 


6.6 


5.4 


10'^' 


5.4 


10'^* 


617 


544 


A1B3G3-D3M11 


58.5 


29.1 


0.68 


8.3 


2.2 


• 10'* 


2.2 


10'* 


330 


2.77 


A1B3G3-D3M12 


70 


23.6 


0.70 


4.6 


1.4 


10'^ 


1.2 


10'5 


625 


557 


A1B3G3-D3M13 


70 


27.1 


0.73 


140.7 


9.9 


10'* 


1.2 


10'* 


1068 


568 


A1B3G5-D3M10 


48 


13.6 


0.21 


4.9 


4.4 


10'^' 


4.4 


10'^ 


282 


278 


A1B3G5-D3M11 


48 


13.6 


0.21 


4.8 


4.9 


10'* 


4.9 


10'* 


282 


278 


A1B3G5-D3M12 


48 


14.2 


0.22 


4.9 


8.8 


10'* 


7.2 


10'* 


282 


278 


A1B3G5-D3M13 


48 


13.6 


0.24 


-215 


1.9 


10'5 


4.8 


10'3 


285 


278 


A2B4G1-D3M10 


114 


119 


0.86 


58.1 


3.9 


•10" 


3.9 


10" 


518 


438 




151 


110 


1.24 


29.4 


3.6 


•10" 


3.6 


10" 








179 


118 


0.91 


57.0 


6.9 


10" 


6.9 


10" 






A2B4G1-D3M11 


115 


119 


0.85 


58.0 


4.0 


10'2 


3.9 


10'2 


518 


438 
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111 


1.24 


29.4 


3.6 


10'2 


3.6 


10'^ 






A2B4G1-D3M12 


124 
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0.76 


72.2 


2.6 


10" 


2.1 


10" 


881 


634 


A2B4G1-D3M13 


123 


120 


1.1 


40.3 


1.5 


10'* 


1.1 


10'* 


570 


506 




143 


57.4 


1.1 


8.3 


5.4 


10'* 


3.3 


10'* 
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60.6 


1.1 


8.2 


5.8 


10" 


3.1 


10" 






A2B4G4-D3M10 


48 


28 


0.47 


6.4 


8.8 


•10'^ 


8.6 


10'^ 


246 


186 




60 


29.5 


0.48 


7.1 


1.1 


•10'3 


1.1 


10'3 


482 


370 


A2B4G4-D3M11 


48 


28 


0.47 


6.4 


8.7 


•10'3 


8.6 


10" 


246 


184 




60 


29.4 


0.48 


7.0 


9.3 


•10'3 


9.0 


10" 


482 


370 


A2B4G4-D3M12 


48 


28.3 


0.47 


6.0 


8.0 


•10'* 


7.8 


10'* 


266 


201 


A2B4G4-D3M13 


48 


17.6 


0.48 


10.6 


1.9 


■ 10'5 


4.5 


10'* 


332 


249 


A2B4G5-D3M10 


50 


13.9 


0.21 


3.1 


2.7 


10'^ 


2.7 


10'^' 


308 


290 


A2B4G5-D3M11 


50 


13.9 


0.22 


3.1 


3.4 


10'* 


3.4 


10'* 


308 


290 


A2B4G5-D3M12 


50 


12.0 


0.23 


3.2 


3.2 


10" 


2.9 


10" 


316 


293 


A2B4G5-D3M13 


50 


15.2 


0.28 


-70.7 


2.4 


10" 


1.1 


10'* 


891 


320 
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'■c 


Mc 


2n/n 


\b\ 


\h\ 




'4 




[ms] 


[km] 


[Mo] 


[ms] 


[G] 


[G] 


[km] 


[km] 


A3B2G4-D3M10 


57 


19.6 


0.44 


4.6 


4.7 


• 10'-' 


4.7 • 


10'' 


378 


340 


A3B2G4-D3M11 


57 


20.0 


0.44 


4.8 


4.0 


• 10'" 


4.0 • 


10'" 


382 


340 


A3B2G4-D3M12 


57 


28.2 


0.44 


8.7 


8.7 


• 10'" 


8.4 • 


10'" 


382 


348 


A3B2G4-D3M13 


57 


17.1 


0.45 


41.2 


2.4 


•10'^ 


2.1 • 


10'" 


906 


405 


A3B3G3-D3M10 


64.7 


54.9 


0.73 


20.3 


2.0 


• 10'^ 


1.9 • 


10'^ 


491 


398 




67.8 


37.4 


0.71 


10.0 


1.1 


• 10'^ 


1.1 • 


10'^ 


564 


463 


A3B3G3-D3M11 


64.7 


54.9 


0.73 


20.3 


2.0 


• 10'^ 


1.9 • 


10'^ 


485 


398 




67.7 


37.5 


0.71 


9.9 


1.2 


■ 10'" 


1.2- 


10'" 


558 


458 


A3B3G3-D3M12 


68 


39.7 


0.65 


6.2 


7.0 


■ 10'" 


6.9 • 


10'" 


564 


469 


A3B3G3-DOM13 


64.2 


63.1 


0.68 


12.8 


4.2 


• 10'" 


2.5 • 


10'" 


911.5 


508.4 


A3B3G3-D1M13 


68 


39.7 


0.71 


14.8 


4.9 


■ 10'3 


4.8 • 


10'3 


571 


469 


A3B3G3-D3M13 


68 


24.7 


0.63 


11.1 


1.1 


• 10'^ 


3.4 • 


10'" 


901 


571 


A3B3G4-D3M10 


61.7 


25.2 


0.47 


6.3 


2.7 


• 10'' 


2.7 • 


10'* 


497 


412 


A3B3G4-D3M11 


62 


24.4 


0.46 


5.7 


4.3 


■ 10'" 


4.3 • 


10'" 


503 


417 


A3B3G4-D3M12 


62 


35.9 


0.44 


5.5 


9.4 


■ 10'" 


7.0 • 


10'" 


520 


417 


A3B3G4-D3M13 


43 


27.72 


0.48 


34.2 


2.3 


• 10^" 


1.4 • 


10'" 


131 


111 


A3B3G5-D3M10 


44.3 


13.4 


0.22 


3.3 


2.8 


• 10'^ 


2.7 • 


10'^* 


224 


221 


A3B3G5-D3M11 


44.3 


13.7 


0.22 


3.6 


3.3 


■ 10'" 


3.3 • 


10'" 


226 


224 


A3B3G5-D0M12 


52.7 


29.2 


0.23 


7.7 


5.3 


• 10'" 


5.2 • 


10'" 


375 


358 


A3B3G5-D1M12 


54.6 


16.4 


0.24 


5.6 


1.3 


• 10'" 


1.3- 


10'" 


371 


371 


A3B3G5-D2M12 


54.6 


13.7 


0.25 


3.4 


2.0 


• 10'5 


1.9- 


10'5 


384 


388 


A3B3G5-D3M12 


56.9 


12.6 


0.24 


2.3 


3.1 


• 10'^ 


2.8- 


10'5 


417 


417 


A3B3G5-D4M12 


57.6 


14.8 


0.24 


3.4 


5.3 


• 10'" 


5.2- 


10'" 


375 


358 


A3B3G5-D3M13 


59.6 


15.1 


0.27 


-38.0 


2.3 


■ 10'5 


8.9- 


10'3 


1009 


354 


A4B5G5-D3M10 


52.7 


86.0 


0.74 


21.6 


2.3 


■10'^ 


1.2- 


10'^ 






A4B5G5-D3M11 


52.7 


71.3 


0.68 


17.1 


2.3 


• 10'3 


2.22 


•10'3 






A4B5G5-D3M12 


36.1 


74.4 


0.70 


10.9 


1.4 


• 10'" 


1.0- 


10'" 






A4B5G5-D3M13 


46.8 


13.9 


0.36 


34.8 


3.4 


• 10'^ 


2.0- 


10'" 







